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Dynamical tunneling, introduced in the molecular context, is more than two decades old and refers 
to phenomena that are classically forbidden but allowed by quantum mechanics. The barriers for 
dynamical tunneling, however, can arise in the momentum or more generally in the full phase space of 
the system. On the other hand the phenomenon of intramolecular vibrational energy redistribution 
(IVR) has occupied a central place in the field of chemical physics for a much longer period of time. 
Despite significant progress in understanding IVR a priori prediction of the pathways and rates is 
still a difficult task. Although the two phenomena seem to be unrelated several studies indicate that 
dynamical tunneling, in terms of its mechanism and timescales, can have important implications 
for IVR. It is natural to associate dynamical tunneling with a purely quantum mechanism of IVR. 
Examples include the observation of local mode doublets, clustering of rotational energy levels, and 
extremely narrow vibrational features in high resolution molecular spectra. Many researchers have 
demonstrated the usefulness of a phase space perspective towards understanding the mechanism of 
IVR. Interestingly dynamical tunneling is also strongly influenced by the nature of the underlying 
classical phase space. Recent studies show that chaos and nonlinear resonances in the phase space 
can enhance or suppress dynamical tunneling by many orders of magnitude. Is it then possible 
that both the classical and quantum mechanisms of IVR, and the potential competition between 
them, can be understood within the phase space perspective? This review focuses on addressing 
the question by providing the current state of understanding of dynamical tunneling from the phase 
space perspective and the consequences for intramolecular vibrational energy flow in polyatomic 
molecules. 
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INTRODUCTION 

The field of chemical physics is full of phenom- 
ena that occur quantum mechanically, with observ- 



able consequences, despite being forbidden by classical 
mechanics [lj]. All such processes are labeled as tunnel- 
ing with the standard elementary example being that of 
a particle surmounting a potential barrier despite insuf- 
ficient energy. However the notion of tunneling is far 
more general in the sense that the barriers can arise in 
the phase space @, EJUS]. In other words the barri- 
ers are dynamical and arise due to the existence of one 
or several conserved quantities. Barriers due to exactly 
conserved quantities i.e., constants of the motion are usu- 
ally easy to identify. A special case is that of a particle 
in a one dimensional double well potential wherein the 
barrier is purely due to the conserved energy. On the 
other hand it is possible, and frequently observed, that 
the dynamics of the system can result in one or more 
approximate constants of the motion which can manifest 
as barriers. The term approximate refers to the fact that 
the relevant quantities, although strictly nonconserved, 
are constant over timescales that are long compared to 
certain system timescales of interest @. Such approxi- 
mate dynamical barriers are not easy to identify and in 
combination with the other exact constants of the mo- 
tion can give rise to fairly complex dynamics. As usual 
one anticipates that the dynamical barriers will act as 
bottlenecks for the classical dynamics whereas quantum 
dynamics will 'break free' due to tunneling through the 
dynamical barriers i.e., dynamical tunneling. However, 
as emphasized in this review, the situation is not nec- 
essarily that straightforward since the mere existence of 
a finite dynamical barrier does not guarantee that dy- 
namical tunneling will occur. This is especially true in 
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multidimensions since a variety of other dynamical ef- 
fects can localize the quantum dynamics. In this sense 
the mechanism of dynamical tunneling is far more sub- 
tle as compared to the mechanism of tunneling through 
nondynamical barriers. 

The distinction between energetic and dynamical 
barriers can be illustrated by considering the two- 
dimensional Hamiltonian 



H(x,y,p x ,p y ) 
V(x,y) 



(Pl+P 2 y ) + V(x,y) 



(1) 



-ir i 2 2\2 i 1 2 2 i 2 2 

V (x - a ) + - uj y + jxy 



which represents a one-dimensional double well (in the x 
degree of freedom) coupled to a harmonic oscillator (in 
the y degree of freedom). For energies E below the bar- 
rier height Vb = Vna 4 , the y — Poincare surface of sec- 
tion in Fig.[T]shows two disconnected regions in the phase 
space. This reflects the fact that the corresponding isoen- 
ergetic surfaces are also disconnected. Thus one speaks 
of an energetic barrier separating the left and the right 
wells. On the other hand Fig.Q]shows that for E > V> the 
Poincare surface of section again exhibits two regular re- 
gions related to the motion in the left and right wells de- 
spite the absence of the energetic barrier. In other words 
the two regions are now part of the same singly connected 
energy surface but the left and the right regular regions 
are separated by dynamical barriers. The various non- 
linear resonances and the chaotic region seen in Fig. Q] 
should be contrasted with the phase space of the one- 
dimensional double well alone which is integrable. Later 
in this review examples will be shown wherein there are 
only dynamical barriers and no static potential barriers. 
Note that in higher dimensions, see discussions towards 
the end of this secion, one does not even have the lux- 
ury of visualizing the global phase space, as in Fig. [TJ let 
alone identifying the dynamical barriers! 

In the case of multidimensional tunneling through po- 
tential barriers it is now well established that valuable in- 
sights can be gained from the phase space perspective @, 
S 0, EE El It is not necessary that tunneling is as- 
sociated with transfer of an atom or group of atoms 
from one site to another site[12]. One can have, for in- 
stance, vibrational excitations transferring from a par- 
ticular mode of a molecule to a completely different 
mode[I3,[ll. Examples like Hydrogen atom transfer and 
electron transfer belong to the former class whereas the 
phenomenon of intramolecular vibrational energy redis- 
tribution (IVR) is associated with the latter class fl3j|. 
In recent years dynamical tunneling has been realized in 
a number of physical systems. Examples include driven 
atoms 14, 15 1, microwave llil or optical cavities [H|. Bose- 
Einstein condensates (lH T2Q] . and in quantum d ots[2lll . 



Thinking of dynamical tunneling as a close cousin IE 
of the "above barrier reflection" (cf. example shown in 
Fig. [T]) a recent paper by Giese et al. suggests the im- 




Figure 1: Poincare surface of section for the example two- 
dimensional potential in Eq. {T}. The parameter values are 
chosen to be Vb = 0.25, a = uj — 7 = 1 and the static barier 
height is therefore V& = 0.25. The surface of section are de- 
fined by y = and p y > 0. The top panel shows the phase 
space for E — 0.2 (below the static barrier height) and one ob- 
serves two disconnected regular regions separated by energetic 
barrier. The bottom panel corresponds to E — 0.4 (above 
the static barrier height) and shows the left and right regu- 
lar regions separated now only by dyn amica l barriers. Model 
system and parameters taken from ref. 1 1791 . 



portance of dynamical tunneling in order to understand 
the eigenstates of the dichlorotropolone molecule. This 
review is concerned with the molecular manifestations of 
dynamical tunneling and specifically on the relevance of 
dynamical tunneling to IVR and the corresponding signa- 
tures in molecular spectra. One of the aims is to reveal 
the detailed phase space description of dynamical tun- 
neling. This, in turn, leads to the identification of the 
key structures in phase space and a universal mechanism 
for dynamical tunneling in all of the systems mentioned 
above. 

In IVR studies, which is of paramount interest to 
chemical dynamics, one is interested in the fate of an 
initial nonstationary excitation in terms of timescales, 
pathways and destinations [H 0, HE HE HI- Will 



the initial energy "hot-spot" redistribute throughout the 
molecule statistically? Alternatively, to borrow a term 
from a recent review [28j], is the observed statisticality 
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only "skin-deep"? The former viewpoint is at the heart 
of one of the most useful theories for reaction rates - the 
RRKM (Rice-Rampsperger-Kassel- Marcus) theory [29j]. 
On the other hand, recent studies [30l] seem to be lean- 
ing more towards the latter viewpoint. As an aside it 
is worth mentioning that IVR in molecules is essentially 
the FPU 31, 12] (Fermi-Pasta-Ulam) problem; only now 



one needs to worry about a multidimensional network of 
coupled nonlinear oscillators. In a broad sense, the hope 
is that a mechanistic understanding of IVR will yield 
important insights into mode-specific chemistry and the 
coherent control of reactions. Consequently, substantial 
experimental 25, 33, 34| and theoretical efforts [35L I3H 
have been directed towards understanding IVR in both 
time and frequency domains. Most of the studies, span- 
ning many decades, have focused on the gas phase. More 
recently, researchers have studied IVR in the condensed 
phase and it appears that the gas phase studies provide a 
useful starting point. A detailed introduction to the liter- 
ature on IVR is beyond the scope of the current article. 
A brief description is provided in the next section and 
the reader is refered to the literature below for a compre- 
hensive account of the recent advances. The review [33j 
by Nesbitt and Field gives an excellent introduction to 
the literature. The review[3|| by Gruebele highlights the 
recent advances and the possibility of controlling IVR. 
The topic of molecular energy flow in solutions has been 
recently reviewed (39| by Afimann, Kling and Abel. 

Tunneling leads to quantum mixing between states lo- 
calized in classically disconnected regions of the phase 
space. In this general setting barriers arise due to ex- 
actly or even approximately conserved quantities. Thus, 
for example, it is possible for two or more localized quan- 
tum states to mix with each other despite the absence of 
any energetic barriers separating them. This has signifi- 
cant consequences for IVR in isolated molecules [13j since 
energy can flow from an initially excited mode to other, 
qualitatively different, modes; classically one would pre- 
dict very little to no energy flow. Hence it is appropriate 
to associate dynamical tunneling with a purely quantum 
mechanism for energy flow in molecules. In order to have 
detailed insights into IVR pathways and rates it is nec- 
essary to study both the classical and quantum routes. 
The division is artificial since both mechanisms coexist 
and compete with each other. However, from the stand- 
point of control of IVR deciding the importance of one 
route over the other can be useful. In molecular sys- 
tems the dynamical barriers to IVR are related to the 
existence of the so called polyad numbers mEl. Usu- 
ally the polyad numbers are quasiconserved quantities 
and act as bottlenecks for energy flow between different 
polyads. Dynamical tunneling, on the otherhand, could 
lead to interpolyad energy flow. Historically the time 
scales for IVR via dynamical tunneling have been thought 
to be of the order of hundreds of picoseconds. However 
recent advances in our understanding of dynamical tun- 



neling suggest that the timescales could be much smaller 
in the mixed phase space regimes. This is mainly due to 
the observations that chaos in the underlying phase 
space can enhance the tunneling by several orders of mag- 
nitude. Conversely, very early on it was thought that the 
extremely long timescales would allow for mode-specific 
chemistry. Presumably the chaotic enhancement would 
spoil the mode specificty and hence render the initially 
prepared state unstable. Thus it is important to under- 
stand the effect of chaos on dynamical tunneling in order 
to suppress the enhancement. 

Obtaining detailed insights into the phenomenon of 
dynamical tunneling and the resulting consequences for 
IVR requires us to address several questions. How does 
one identify the barriers? Can one explicitly construct 
the dynamical barriers for a given system? What is the 
role of the various phase space structures like resonances, 
partial barriers due to broken separatrices and cantori, 
and chaos? Which, if any, of the phase space structures 
provide for a universal description and mechanism of dy- 
namical tunneling? Finally the most important question: 
Do experiments encode the sensitivity of dynamical tun- 
neling to the nature of the phase space? This review is 
concerned with addressing most of the questions posed 
above in a molecular context. Due to the nature of the 
issues involved considerable work has been done by the 
nonlinear physics community and in this work some of the 
recent developements will be highlighted since they can 
have significant impact on the molecular systems as well. 
Another reason for bringing together developements in 
different fields has to do with the observation that the 
literature on dynamical tunneling and its consequences 
are, curiuosly enough, disjoint with hardly any overlap 
between the molecular and the 'non' molecular areas. 

At this stage it is important and appropriate, given the 
viewpoint adopted in this review, to highlight certain cru- 
cial issues pertaining to the nature of the classical phase 
space for systems with several (N > 3) degrees of free- 
dom. A significant portion of the present review is con- 
cerned with the phase space viewpoint of dynamical tun- 
neling in systems with N < 2. The last section of the re- 
view discusses recent work on a model with three degrees 
of freedom. The sparsity of work in N > 3 is not entirely 
surprising and parallels the situation that prevails in the 
classical-quantum correspondence studies of IVR in poly- 
atomic molecules. Indeed it will become clear from this 
review that classical phase space structures are, to a large 
extent, responsible for both classical and quantum mech- 
anisms of IVR. Thus, from a fundamental, and certainly 
from the molecular, standpoint it is highly desirable to 
understand the mechanism of classical phase space trans- 
port in systems with N > 3. In the early eighties the sem- 
inal work [111 by Mackay, Meiss, and Percival on trans- 
port in Hamiltonian systems with N = 2 motivated re- 
searchers in the chemical physics community to investi- 
gate the role of various phase space structures in IVR 
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dynamics. In particular these studies 44j, El H Eil Ej 
helped in providing a deeper understanding of the dy- 
namical origin of nonstatistical behaviour in molecules. 
However molecular systems typically have atleast N = 3 
and soon the need for a generalization to higher degrees 
of freedom was felt; this task was quite difficult since 
the concept of transport across broken separatrices and 
cantori do not have a straightforward generalization in 
higher degrees of freedom. In addition tools like the 
Poincare surface of section cease to be useful for visual- 
izing the global phase space structures. Wiggins [49| pro- 
vided one of the early generalizations based on the con- 
cept of normally hyperbolic invariant manifolds (NHIM) 
which led to illuminating studies I5H in order to elu- 
cidate the dynamical nature of the intramolecular bot- 
tlenecks for N > 3. Although useful insights were gained 



from several studics|50, 51, 52, \5& 54, 55, \m the 



tramolecular bottlenecks could not be characterized at 
the same levels of detail as in the two degree of free- 
dom cases. It is not feasible to review the various studies 
and their relation/consequences to the topic of this ar- 
ticle. The reader is referred to the paper [50| by Gillilan 
and Ezra for an introduction and to the monograph [57jj 
by Wiggins for an exposition to the theory of NHIMs. 
Following the initial studies, which were perhaps ahead 
of their time, far fewer efforts were made for almost a 
decade but there has been a renewal of interest in the 
problem over the last few years with the NHIMs play- 
ing a crucial role. Armed with the understanding that 
the notion of partial barriers, chaos, and resonances are 
very different in N > 3 fresh insights on transition state 
and RRKM theories, and hence IVR, are beginning to 
emerged! |H 0, M, H, M, 0, M, M, H3, H& HI . 

To put the issues raised above in perspective for the 
current review note that there are barriers in phase space 
through which dynamical tunneling occurs and at the 
same time there are barriers that can also lead to the 
localization of the quantum dynamics. The competition 
between dynamical tunneling and dynamical localization 
is already important for systems with N ~ 2 and this will 
be briefly discussed in the later sections. Towards the end 
of the review some work on a three degrees of freedom 
model are presented. However, as mentioned above, the 
ideas are necessarily preliminary since to date one does 
not have a detailed understanding of either the classical 
barriers to transport nor the dynamical barriers which 
result in dynamical tunneling. Naturally, the competi- 
tion between them is an open problem. We begin with 
a brief overview of IVR and the explicit connections to 
dynamical tunneling. 



State space model of IVR 

Consider a molecule with N atoms which has s = 
3N — 6 (3N — 5 for linear molecules) vibrational modes. 



Dynamical studies, classical and/or quantum, require 
the s-dimensional Born-Oppenheimer potential energy 
surface V(Qi, Q2, ■ ■ ■ , Qs) in terms of some convenient 
generalized coordinates Q and their choice is a cru- 
cial issue. Obtaining the global V(Q) by solving the 
Schrodingcr equation is a formidable task even for mid- 
sized molecules. Traditionally, therefore, a perturbative 
viewpoint is adopted which has enjoyed a considerable 
degree of success. For example, near the bottom of the 
well and for low levels of excitations the approximation of 
vibrations by uncoupled harmonic normal modes is suffi- 
cient and the molecular vibrational Hamiltonian can be 
expressed perturbatively as 7^, 71 | : 



H 



i=i 



Ql) 



^ 4>ijkQiQjQk 
ijk—1 



(2) 



In the above (P, Q) are the dimensionless vibrational mo- 
menta and normal mode coordinates respectively. The 
deviations from the harmonic limit are captured by the 
anharmonic terms with strengths {4>ijk}- Such small 
anharmonicities account for the relatively weak over- 
tone and combination transitions observed in experimen- 
tal spectra. However with increasing energy the anhar- 
monic terms become important and doubts arise regard- 
ing the appropriateness of a perturbative approach. The 
low energy normal modes get coupled and the very con- 
cept of a mode becomes ambiguous. There exists suffi- 
cient evidence 72, 73, 74, 75| for the appearance of new 
'modes', unrelated to the normal modes, with increasing 
vibrational excitations. Nevertheless detailed theoretical 
studies over the last couple of decades has shown that it is 
still possible to understand the vibrational dynamics via 
a suitably generalized perturbative approach - the canon- 
ical Van-Vleck perturbation theory (CVPT) 0, 0, [H . 
The CVPT leads to an effective or spectroscopic Hamil- 
tonian which can be written down as 

H = H Q (n) + V(k,eJ) (3) 

where n = at a, a, and aJ are the harmonic number, de- 
struction and creation operators. The zeroth-order part 
i.e., the Dunham expansion [79j 



H (n) = UjNj + Y x lJ N l N j + 



(4) 



>>j 



is diagonal in the number representation. In the above 
expression Nj = hj + dj/2 with rij being the number 
of quanta in the j th mode with degeneracy dj. The off- 
diagonal terms i.e., anharmonic resonances have the form 



m j 



with m = {mf,'i 



(5) 

, .}. These terms represent the 
couplings between the zeroth-order modes and are re- 
sponsible for IVR starting from a general initial state 
l*)=E J c nj |n). 
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Figure 2: Schematic diagram illustrating the tier model of 
IVR. The bright state |6) is coupled to the various dark states 
sorted into tiers through perturbations V. As an example 
this sketch indicates fast IVR between |6) and states in |d)i 
whereas further flow of energy to states \d)2 in the second tier 
is restricted due to the lack of couplings. Nevertheless it is 
possible that vibrational superexchange can couple the first 
and second tiers and energy can eventually flow into the final 
tier (grey) over long timescales. 



A few key points regarding the effective Hamiltonians 
can be noted at this juncture. There are two routes to 
the effective Hamiltonians. In the CVPT method the 
parameters of the effective Hamiltonian are related to 
the original molecular parameters of V(Q). The other 
route, given the difficulties associated with determining 
a sufficiently accurate I^(Q), is via high resolution spec- 
troscopy. The parameters of Eq.Q are determined (70| 
by fitting the experimental data on line positions and, to 
a lesser extent, intensities. The resulting effective Hamil- 
tonian is only a model and hence the parameters of the 
effective Hamiltonian do not have any obvious relation- 
ship to the molecular parameters of V(Q). Although the 
two routes to Eq.((3]) are very different, they complement 
one another and each has its own advantages and draw- 
backs. The reader is refered to the reviews by Sibert and 
Joyeux for a detailed discussion of the CVPT approach. 
The review 72] by Ishikawa et al. on the experimental 
and theoretical studies of the HCP <-*■ CPH isomeriza- 
tion dynamics provides, amongst other things, a detailed 
comparison between the two routes. 

Imagine preparing a specific initial state, say an eigen- 
state of Ho, denoted as |b). Theoretically one is 
not restricted to |b) and very general class of initial 
states can be considered. However, time domain exper- 
iments with short laser pulses excite specific vibrational 
overtone or combination states, called as the zeroth- 
order bright states (ZOBS), which approximately corre- 
spond to the eigenstates of Hq. The rest of the opti- 



cally inaccessible states are called as dark states {|d)}. 
Here the zeroth-order eigenstates |n) are partitioned as 
{|b), |d)i, |d) 2 , . . .}. The perturbations V(a, a*) couple 
the bright state with the dark states leading to energy 
flow and impose a hierarchical coupling structure be- 
tween the ZOBS and the dark states. Thus one imagines, 
as shown in Fig.[2j the dark states to be arranged in tiers, 



determined by the order p = J2i( m j 



m. 



-) of the V a 



with an increasing density of states across the tiers. This 
implies that the the local density of states around is 
the important factor and experiments do point to a hi- 
erarchical IVR process 23j. For example, Callegari et al. 
have observed seven different time scales ranging from 
100 femtoseconds to 2 nanoseconds for IVR out of the 
first CH stretch overtone of the benzene molecule [80jj . 
Grucbele and coworkers have impressively modeled the 
IVR in many large molecules based on the hierarchical 
tiers concept Q. In a way the coupling V(a, a*) imposes 
a tier structure on the IVR; similar tiered IVR flow would 
have been observed if one had investigated the direct dy- 
namics ensuing from the global V(Q). Thus the CVPT 
and its classical analog help in unraveling the tier struc- 
ture. Nevertheless dominant classical and quantum en- 
ergy flow routes still have to be identified for detailed 
mechanistic insights. 

In the time domain the survival probability 



P h (t) = \(b(0)\b(t))f 



a,/3 



(6) 



gives important information on the IVR process. The 
eigenstates of the full Hamiltonian have been denoted 
by | a), |/3), . . . with pb,a = |(a|b)| 2 being the spectral 
intensities. The long time average of Pb(t) 



<r b = lim I f P b (t)dt = 5>b,, 



(7) 



is known as the inverse participation ratio (also called 
as the dilution factor in the spectroscopic community). 
Essentially a^ 1 indicates the number of states that par- 
ticipate in the IVR dynamics out of |b). 

Although the tier picture has proved to be useful in 
analysing IVR, recent studies 25, 8l|, [H, HH, HI] empha- 
size the far greater utility in visualising the IVR dynam- 
ics of |b(t)) as a diffusion in the zeroth-order quantum 
number space or state space denoted by (ni, n^, ■ ■ ■ , n s ). 
The state space picture, shown schematically in Fig. [3j 
highlights the local nature and the directionality of the 
energy flow due to the various anharmonic resonances. 
In the "one-dimensional" tier scheme shown in Fig. [5] 
the anisotropic nature of the IVR dynamics is hard to 
discern. Tiers in the state space are formed by zeroth- 
order states within a certain distance from the bright 
state and hence still organized by the order of the cou- 
pling resonances. In the long time limit most states that 
participate in the IVR dynamics are confined, due to the 
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Leitner and Wolynes^ building upon the earlier work 
of Logan and Wolynes[86|, have provided criteria[87| for 
vibrational state mixing and energy flow from the state 
space perspective. Using a local random matrix approach 
to the Hamiltonian in Eq. ([3]) the rate of energy flow out 
of |b) is given by 



(9) 



with Q being a distance in the state space. The term 
(|V>q|) represents the average effective coupling of |b) to 
the states, Kq of them with density Dq(E), a distance 
Q away in state space. The extent of state mixing is 
characterized by the transition parameter 



Figure 3: (Colour online) Schematic diagram illustrating the 
state space model of IVR. Zeroth order states on the energy 
shell (blue circles) and off the energy shell (green circles) are 
indicated. An interior state (large red circle) can couple to 
nearby states via several anharmonic resonances indicated by 
purple arrows. Edge states (shaded red circle) can couple via 
off-resonant virtual states \i a ) to states on the energy shell. 
Examples of such one-step and multistep superexchange paths 
in the state space are shown. 



T W = Y \J2 K Q(^Q\) D Q( E )j (10) 

and the transition between localized and extended states 
is located at T(E) = 1. The key term is the effective 
coupling ipQ which involves both low and high order res- 
onances. Applications to several systems shows that the 
predictions based on Eq. ^ and Eq. (flT)]) are both qual- 
itatively and quantitatively [88{ accurate. 



quasi-microcanonical nature of the laser excitation, to 
an energy shell E® ± SE with E® being the bright state 
energy. The IVR dynamics occurs on a _D-dimcnsional 
subspace, refered to as the IVR "manifold" , of the state 
space due to the local nature of the couplings between |b) 
and |dj) with H i=s E® w E%.. One of the key prediction 
of the state space model is that the survival probability 
exhibits a power law at intermediate times (271. |85| 



Pb(t) ~ <r b + (1 - <7b) 



1 



2t 



-D/2 



(8) 



with D <C 37V— 6, the number of vibrational modes of the 
molecule. Thus the effective dimensionality D of the IVR 
manifold in the state space can be fairly small and a large 
body of work in the recent years support the state space 
viewpoint 2f| 27, 3j|. The effective dimensionality itself 
is crucially dependent on the extent to which classical 
and quantum mechanisms of IVR manifest themselves at 
specific energies. One possible interpretation of D is as 
follows. If D 3N — 6 then the dynamics in the state 
space can be thought of as a normal diffusive process and 
thus ergodic over the state space. In such a limit for large 
N (molecules) the survival can be well approximated by 
an exponential behaviour. For D <C 3 V — 6 the IVR 
dynamics is anisotropic and the dynamics is nonergodic 
with D typically being nonintegral. Note that the terms 
"manifold" and "effective dimensionality" are being used 
a bit loosely since one does not have a very clear idea of 
the topology of the IVR manifold as of yet. 



Connections to dynamical tunneling 

The perspective of IVR being a random walk on an ef- 
fective D-dimensional manifold in the quantum number 
space is reasonable 25, 81] as long as direct anharmonic 
resonances exist which connect the bright state with the 
dark states. It is useful to emphasize again that the ex- 
istence of direct anharmonic resonances in itself does not 
imply ergodicity over state space i.e., the random walk 
need not be normal. What happens if, for certain bright 
states, there are no direct resonances available? In other 
words, the coupling matrix elements (dj|V|b) w for 
all j. What would be the mechanism of IVR, if any, in 
such cases? Examples of such states in fact correspond 
to overtone excitations which are typically prepared in 
experiments. The overtone states are also called as edge 
states from the state space viewpoint since most of the 
excitation is localized in one single mode of the molecule. 
Thus, for example, in a system with four degrees of free- 
dom a state |b) = |m, 0, 0, 0) would be called as an edge 
state whereas the state |m, 0, 713, 714), corresponding to a 
combination state, is called as an interior state. The edge 
states, owing to their location in the state space, have 
fewer anharmonic resonances available for IVR as com- 
pared to the interior states. To some extent such a line 
of reasoning leads one to anticipate overtone excitations 
to have slower IVR when compared to the combination 
states. A concrete example comes from the theoretical 
work of Holme and Hutchinson wherein the dynamics 
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of overtone excitations in a model system representing 
coupled CH-stretch and CCH-bend interacting with an 
intense laser field was studied [H, They found that 
classical dynamics predicted no significant energy flow 
from the high overtone excitations. However the corre- 
sponding quantum calculations did indicate significant 
IVR. Based on their studies Holme and Hutchinson con- 
cluded that overtone absorption proceeds via dynamical 
tunneling on timescales of the order of a few nanoseconds. 

Experimental evidence for the dynamical tunneling 
route to IVR was provided in a series of elegant works [9 ll 
|92| by the Princeton group. The frequency domain ex- 
periments involved the measurement of IVR lifetimes of 
the acetylinic CH-stretching states in (CX3)3YCC=H 
molecules with X=H,D and Y=C,Si. The homogeneous 
line widths, which are related to the rate of IVR out of 
the initial nonstationary state, arise due to the vibra- 
tional couplings between the CH-stretch and the various 
other vibrational modes of the molecule. Surprisingly 
Kerstel et al. found [9l| extremely narrow linewidths 



of the order of 10 



10 cm which translate to 



timescales of the order of thousands of vibrational pe- 
riods of the CH-stretching vibration. Thus the initially 
prepared CH-stretch excitation remains localized for ex- 
tremely long times. Such long IVR timescales were as- 
cribed to the lack of strong/direct anharmonic resonances 
coupling the CH-stretch with the rest of the molecular 
vibrations. The lack of resonances, despite a substan- 
tial density of states, combined with IVR timescales of 
several nanoseconds implies that the mechanism of en- 
ergy flow is inherently quantum. Another studyj^ by 
Gambogi et al. examines the possibility of long range 
resonant vibrational energy exchange between equivalent 
CH-stretches in CH3Si(C=CH)3. There have been sev- 
eral experiments [93, [H, @, 97, [98| that indicate mul- 
tiquantum zeroth-order state mixings due to extremely 
weak couplings of the order oflO -1 - 10~ 2 cm -1 in highly 
vibrationally excited molecules. 

Since the early work by Kerstel et al. other exper- 
imental studies have revealed the existence of the dy- 
namical tunneling mechanism for IVR in large organic 
molecules. For instance in a recent work[99] Callegari 
et al. performed experimental and theoretical studies 
on the IVR dynamics in pyrrole (C4H4NH) and 1,2,3- 
triazine (C3H3N3). Specifically, they chose the initial 
bright states to be the edge state 2i^i 4 CH-stretch for 
pyrrole and the interior state v\ + corresponding 
to the CH stretching-ring breathing combination for tri- 
azine. In both cases very narrow IVR features, similar 
to the observations by Kerstel et al., were seen in the 
spectrum. This pointed towards an important role of the 
off-resonant states to the observed and calculated narrow 
IVR features. The analysis by Callegari et al. reveals 
that near the IVR threshold it is reasonable to expect 
such highly off-resonant coupling mechanisms to be op- 
erative. Another example for the possible existence of the 



off- resonant mechanism comes from the experiment [lOfJ 
by Portonov, Rosenwaks, and Bar wherein the IVR dy- 
namics ensuing from 2v\, 3^i, and Ai>i CH- acetylinic 
stretch of 1-butyne are studied. It was found that the 
homogeneous linewidth of the Ai/i state, ~ 0.5 cm -1 , is 
about a factor of two smaller than the widths of 2vi and 
3ui. 

It is important to note that the involvement of such 
off-resonant states can occur at any stage of the IVR 
process from a tier perspective. In other words it is 
possible that the bright state undergoes fast initial IVR 
due to strong anharmonic resonances with states in the 
first tier but the subsequent energy flow might be ex- 
tremely slow. Several experiments point towards the ex- 
istence of such unusually long secondary IVR time scales 
which might have profound consequences for the interpre- 
tation of the high r esolution spectra. Boyarkin and Rizzo 
demonstrated 101| the slow secondary time scales in their 
experiments on the IVR from the alkyl CH-stretch over- 
tones of CF3H. Upon excitation of the CH-stretch fast 
IVR occurs to the CH stretch-bend combination states 
on femtosecond timescales. However the vibrational en- 
ergy remains localized in the mixed stretch-bend states 
and flows out to the rest of the molecule on timescales of 
the order of hu ndred picoseconds. Similar observations 
were made 102] by Lubich et al. on the IVR dynamics of 
OH-stretch overtone states in CH3OH. 

The interpretation of the above experimental results as 
due to dynamical tunneling is motivated by the theoreti- 
cal analysis by Stuchebrukhov and Marcus in a landmark 
paperS In the initial work they explained the nar- 
row features, observed in the experiments [2H by Kerstel 
et al., by invoking the coupling of the bright states with 
highly off-resonant gateway states which, in turn, cou- 
ple back to states that are nearly isoenergetic with the 
bright state. In Fig. [3] an example of the indirect cou- 
pling via a state which is off the energy shell in the state 
space is illustrated. More importantly Stuchebrukhov 
and Marcus argued that the mechanism, described a little 
later in this review and a nticipated earlier by Hutchinson, 
Sibert, and Hvnes 1041 ] . involving high order co uplin g 
chains is essentially a form of generalized tunneling [105j. 
Hence one imagines dynamical barriers which act to pre- 
vent classical flow of energy out of the initial CH-stretch 
whereas quantum dynamical tunneling does lead to en- 
ergy flow, albeit very slowly. Similar argum ents had been 
put forward by Hutchinson demonstrating [106| the im- 
portance of the off-resonant coupling leading to the mix- 
ing of nearly degenerate high energy zeroth-order states 
in cyanoacetylene HCC=CN. It was argued that such 
purely quantum relaxation pathway would explain the 
observed broadening of the overtone bands in various 
substituted acetylenes. 

The above examples highlight the possible connections 
between IVR and dynamical tunneling. However it was 
realized very early on that observations of local mode 
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doublets [107|, 11081 1103 ] in molecular spectra and the 
clustering of rotational energy sublevels with high an- 
gular momenta^ can also be associated with dynami- 
cal tunneling. Indeed much of our understanding of the 
phenomenon of dynamical tunneling in the context of 
molecular spectra comes from these early works. One 
of the key feature of the analysis was the focus on a 
phase space description of dynamical tunneling i.e., iden- 
tifying the disconnected regions in the phase space and 
hence classical structures which could be identified as 
dynamical barriers. Thus the mechanism of dynamical 
tunneling, confirmed by computing the splittings based 
on specific phase space structures, could be understood in 
exquisite detail. However such detailed phase space anal- 
ysis are exceedingly difficult, if not impossible, for large 
molecules. For example in the case of the (CXa^YCCH 
system there are 3A — 6 = 42 vibrational degrees of free- 
dom. Any attempt to answer the questions on the origin 
and location of the dynamical barriers in the phase space 
is seemingly a futile excercise. It is therefore not sur- 
prising that Stuchebrukhov and Marcus, although aware 
of the phase space perspective, provided a purely quan- 
tum explanation for t he IVR in terms of the chain of 
off- resonant states [l03j ]. 

Thus one cannot help but ask if the phase space view- 
point is really needed. The answer is affirmative and 
many reasons can be provided that justify the need and 
utility of a phase space viewpoint. First, as noted by 
Heller 0, the concept of tunneling is meaningless 
without the classical mechanics of the system as a base- 
line. In other words, in order to label a process as purely 
quantum it is imperative that one establishes the ab- 
sence of the process within classical dynamics. Note that 
for dynamical tunneling it might not be easy to a priori 
make such a distinction. There are mechanisms of clas- 
sical transport, especially in systems with three or more 
degrees of freedom, involving long timescales which might 
give an impression of a tunneling process. It is certainly 
not possible to differentiate between classically allowed 
and forbidden mechanisms by studying the spectral in- 
tensity and splitting patterns alone due to the nontriv- 
ial role played by the stochasticity in the classical phase 
space. Secondly, the quantum explanation is in essence a 
method to calculate the contribution of dynamical tun- 
neling to the IVR rates in polyatomic molecules. In or- 
der to have a complete qualitative picture of dynami- 
cal tunneling it is n ecess ary, as emphasised by Stuche- 
brukhov and Marcus |l05j . to find an explicit form of the 
effective dynamical potential in the state space of the 
molecule. Third reason, mainly due to the important 
insights gained from recent studies, has to do with the 
sensitivity of dynamical tunneling to the stochasticity in 
the phase space [42]. Chaos can enhance as well as su- 
press dynamical tunneling and for large molecules the 
phase space can be mixed regular-chaotic even at ener- 
gies corresponding to the first or second overtone lev- 



els. Undoubtedly the signatures should be present in the 
splitting and intensity patterns in a high resolution spec- 
trum. However the precise mechanism which governs the 
substantial enhancement/suppression of dynamical tun- 
neling, and perhaps IVR rates is not yet clear. Never- 
theless recent developements indicate that a phase space 
approach is capable of providing both qualitative and 
quantitative insights into the problem. 

A final argument in favor of a phase space perspective 
needs to be mentioned. It might appear that the dimen- 
sionality issue is not all that restrictive for the quantum 
studies a la Stuchebrukhov and Marcus. However it will 
become clear from the discussions in the later sections 
that in the case of large systems even calculating, for 
example local mode splittings via the minimal perturba- 
tive expression in Eq. (|56|) can be prohibitively difficult. 
One might argue that a clever basis would reduce the 
number of perturbative terms that need to be considered 
or alternatively one can look for and formulate criteria 
that would allow for a reliable estimate of the splitting. 
Unfortunately, a priori knowledge of the clever basis or 
the optimal set of terms to be retained in Eq. (|56|) im- 
plies a certain level of insight into the dynamics of the 
system. The main theme of this article is to convey the 
message that such insights truly originate from the clas- 
sical phase space viewpoint. With the above reasons in 
mind the next section provides brief reviews of the ear- 
lier approaches to dynamical tunneling from the phase 
space perspective. This will then set the stage to discuss 
the more recent advances and the issues involved in the 
molecular context. 



DYNAMICAL TUNNELING: EARLY WORK 
FROM THE PHASE SPACE PERSPECTIVE 



In a series of pioneering papers 107, lOSi 109l ] . nearly 
a quarter of a century ago, Lawton and Child showed 
that the experimentally observed local mode doublets in 
H2O could be associated with a generalized tunneling in 
the momentum space. In molecular spectroscopy it was 
appreciated from very early on that highly excited spec- 
tra associated with X-H stretching vibrations are bet- 
ter described in terms of l ocal modes ra ther than the 
conventional normal modes [lloL 111 . 112 1. In a zeroth- 



order description this corresponds to very weakly cou- 
pled, if not uncoupled, anharmonic oscillators. Every 
anharmonic oscillator, modelled by an appropriate Morse 
function, represents a specific local stretching mode of 
the molecule. The central question that Lawton and 
Child asked was that to what extent are the molecular 
vibrations actually localized in the individual bonds (lo- 
cal modes)? Analyzing the classical dynamics on the 
Sorbie-Murrell potential energy surface for H2O, focus- 
ing on the two local O-H stretches, it was found that a 
clear distinction between normal mode and local mode 



9 



behaviour could be observed in the phase space. The 
classical phase space at appropriate energies revealed Q 
two equivalent but classically disconnected regions as a 
signature of local mode dynamics. Lawton and Child im- 
mediately noted the topological similarity between their 
two dimensional Poincare surface of section and the phase 
space of a one dimensional symmetric double well poten- 
tial. However they also noted that the barrier was in 
momentum space and hence the lifting of the local mode 
degeneracy was due to a generalized tunneling in th e mo- 
mentum spa ce. D etailed classical [2J, quantum [lOTf, and 
semiclassical |l08l | investigations of the system allowed 
Lawton and Child to provide an approximate formula 
for the splitting AE nm between two symmetry related 
local mode states \nm ± ) as: 



2huj r , 



■ exp 



|P| ■ dQ 



(11) 




The variables (P,Q) corres pond to the mass-weighted 
momenta and coordinates [1, 107 1. It is noteworthy that 
Lawton and Child realized, correctly, that the frequency 
factor Tiuj nm should be evaluated at a fixed total quantum 
number v = n + m. However the choice of the tunnel- 
ing path C proved to be more difficult. Eventually C 
was taken to be a one dimensional, nondynamical path 
in the two dimensional coordinate space of the symmet- 
ric and antisymmetric normal modes. Although Eq. (fTTj) 
correctly predicts the trend of decreasing AE nm for in- 
creasing m and fixed n, the origin and properties of the 
dynamical barrier remained unclear. 

The analysis by Lawton and Child, although specific 
to the vibrational dynamics of H2O, suggested that the 
phenomenon of dynamical tunneling could be more gen- 
eral. This was confirmed in an influential paper by Davis 
and Heller wherein it was argued that dynamical tun- 
neling could have significant effects on bound states of 
polyatomic molecules Q. Incidentally, the word dynam- 
ical tunneling was first introduced by Davis and Heller. 
This work is remarkable in man y resp ects and provided 
a fairly detailed correspondence [ll5| between dynami- 
cal tunneling and the structure of the underlying clas- 
sical phase space for the Hamiltonian H(s,p s ,u,p u ) — 
(pI +P«)/2 + V(s, u). Specifically, the analysis was done 
on the following two-dimensional potential: 



7(*,«) = i(o;> 



+ u) 2 u u 2 ) + Xu 2 



(12) 



The system has a discrete symmetry V(s, —it) — V(s, u). 
For the specific choice of parameter values uj s = 1.0, uj u = 
1.1, A = — 0.11, and K = 1 the dissociation energy is equal 
to 15.125 and the potential supports about 115 bound 
states. At low energies no doublets are found. How- 
ever above a certain energy, despite the lack of an en- 
ergetic barrier, near-degenerate eigenstates 4>i,4>2 were 
observed with small splittings AE. In analogy with 



Figure 4: (Colour online) An example of the local mode dou- 
blet in the Davis-Heller system. In (a) the (s,u) space rep- 
resentation of the even state (E = 13.589216) is shown. In 
(b) and (c) the linear combination ip± = (<j>i ± 4>2)/V2 are 
shown. Right below the configuration space representations 
the corresponding phase space Husimi representations in the 
(u,p u ) Poincare surface of section are shown. The rightmost 
figure on the top panel shows the dynamical tunneling oc- 
curing when starting with the initial nonstationary state ip+ . 
The period of the oscillations is consistent with the splitting 
AE » 9.125 x 10" 5 . The rightmost fi gure in the bottom panel 
shows the (u,p n ) surface of section. 



the symmetric double well system, linear combinations 
ip± = (0i ±</>2)/v2 yielded states localized in the config- 
uration space (s,u). One such example is shown in Fig. [4] 
In the same figure the probabilities \(ip+(0)\ip+(t))\ 2 and 
\(ip-(0)\ip + (t))\ 2 are also shown confirming the two state 
scenario. The main issue once again had to do with 
the origin and nature of the barrier. Davis and Heller 
showed ■ 



115| that the doublets could be associated with 



the formation of two classically disconnected regions in 
the phase space - very similar to the observations by Law- 
ton and ChildQ. The symmetric stretch periodic or- 
bit, stable at low energies, becomes unstable leading to 
the topological change in the phase space; a separatrix 
is created which separates the normal and local mode 
dynamics. In Fig. [5] the phase space are shown for in- 
creasing total energy and the creation of the separatrix 
can be clearly seen in the {u,p u ) surface of section. The 
near-degenerate eigenstates correspond to the local mode 
regions of the phase space and several such pairs were 
identified upto the dissociation energy of the system. The 
corresponding phase space representation shown in Fig. [4] 
(bottom panel) confirms the localized nature of ip± . Thus 
a key structure in the phase space, a separatrix, arising 
due to a 1:1 resonance between the two modes is respon- 
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sible for the dynamical tunneling. Although Davis and 
Heller gave compelling arguments for the connections be- 
tween phase space structures and dynamical tunneling, 
explicit determination of the dynamical barriers and the 
resulting splittings was not attempted. 
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Figure 5: Poincare surface of sections for the Davis-Heller 
Hamiltonian. The top panel shows the (u,p u ) section and 
the bottom panel shows the {s,p„) surface of section for in- 
creasing total energy E = 5.0 (A), 9.0 (B), and 14.0 (C). 
Note that the symmetric stretch periodic orbit (u = p u — 0) 
becomes unstable around E « 6.5 whereas the asymmetric 
stretch periodic orbit (s = p s = 0) becomes unstable around 
E^ 8.5. 



The observation by Davis and Heller, implicily present 
in the work by Lawton and Child, regarding the impor- 
tance of the 1:1 resonance to the dynamical tunneling 
was put to test in an elegant series of papers |l lot Ill7l | 
by Sibert, Rcinhardt, and Hynes. T hese authors inves- 
tigated i n gre at detail the classical [III] and quantum 
dynamics [117} of energy transfer between bonds in ABA 
type triatomics and gave explicit expressions for the split- 
tings. In particular they studied the model H 2 Hamil- 
tonian: 



ff(x,p)= [ \ G 3iPj+U{x S Y\ +G12P1P2 



(13) 



with (2:1,0:2) and (pi,P2) denoting the two OH bond 
stretches (local modes) and their corresponding momenta 
respectively. The bending motion was ignored and the 
stretches were modelled by Morse oscillators 



\rij) are eigenstates of the unperturbed ] th Morse oscilla- 
tor. Due to the coupling the initial states have nontrivial 
dynamics and thus mix with other zeroth-order states. 
An initial state \n%, n 2 ) implies a certain energy distribu- 
tion in the molecule and the nontrivial time dependence 
signals flow of energy in the molecule. Late r studies by 
Hutchinson, Sibert, and Hynes showed 104} that based 
on the dynamics the various initial states could be clas- 
sified as normal modes or local modes. Morever at the 
border between the two classes were initial states that led 
to a 'nonclassical' energy flow. A representative exam- 
ple, using a different but equivalent system, for the three 
classes is shown in Fig. [5] The normal mode regime illus- 
trated with |ni,n 2 ) = |3, 2) shows complete energy flow 
between the modes on the time scale of a few vibrational 
periods - something that occurs classically as well. The 
local mode behaviour for the initial state 1 3, 0) exhibits 
complete energy transfer resulting in the state |0, 3) but 
the timescale is now hundreds of vibrational periods. De- 
spite the large timescales involved it is important to note 
that such a process does not happen classically. Thus 
this is an example of dynamical tunneling. The so called 
nonclassical case illustrated with an initial state |3, 1) is 
also an example of dynamical tunneling. However the 
associated timescale is n owhe re as large when compared 
to the local mode regime |l04l |. 

The vibrational energy transfer process illustrated 
through the initial state |3, 0) and |3, 1) are examples of 
pure quantum routes to en ergy flow. Hutchinson, Sibert, 
and Hynes proposed 104} that the mechanism for this 
quantum energy flow can be understood as an indirect 
state-to-state flow of probability involving normal mode 
intermediate states. For instance, in the case involving 
the initial state |3, 0) the following mechanism was pro- 
posed: 



|3,0)^|2,1} ^|1,2} ^|0,3> 



(15) 



The reason for this indirect route has to do with the 
fact that estimating the splitting directly (at first or- 
der) A = 2(3, OIG12P1P2IO, 3) yields a value which is more 
than an o rder of magnitude smaller than the the actual 
value |l04l |. Indeed note that the indirect route corre- 
sponds to a third order perturbation in the coupling and 
hence it is possible to estimate the contribution to the 
splitting as: 



A (3,0|F|2,l)<2 ) l|y|l,2)(l,2|F|0,3) 



(Elo 



El,? 



(16) 



U{x) = D[l-e- ax X 



(14) 



with D being the OH bond dissociation energy. Due to 
the equivalence of the OH stretches Gn = G22 and the 
coupling strength G12 is much smaller than G\\. The au- 
thors studied the dynamics of initial states \ni, n 2 ) where 



with V = G12P1P2 and E® „ 2 are the zeroth-order ener- 
gies associated with the states |ni,n 2 ). 

Once again a clear interpretation of the mechanism 
comes from ta king a clo ser look at the underlying classi- 
cal phase space [l04 , 116} . Using the classical action-angle 
variables (I, 6) for Morse oscillators it is possible to write 
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the original Hamiltonian in Eq. (fT3|) as: 

n 2 i 2 



H(I,0)= J2 



AD 



V(I,0) (17) 



with the harmonic frequency Q — \j2DG\\a 2 . The 
coupling term V(I, 6) has infinitely many terms of the 
form V pq (T) cos(p9\ — qB-i) which represent nonlinear res- 
onances pu>i ss quj2 involving the nonlinear mode frequen- 
cies LJi = O — n 2 Ii/2D. Arguments were provided [116| 
for the dominance of the p = q = 1 resonance and hence 
to an excellent approximation 



H(I,0) « ]T f 07 ' - 



4D 



fl 2 P 2 



HP 



2^2 



n 2 p 



8D V 813 
H M (P)-H R (p,iP) 



Vii(I) cos(0i-0 2 ) 

V&(P,p) cos(2^) X 
(18) 



The form Hm — Hr originates via a canonical trans- 
formation from (1,0) — > (P, (f>,p,ip) with I\ + I2 — 
P,h ~ h = PA +6 2 = 20, and 6>i - 2 = 2^. Thus 
the key term in the above Hamiltonian is the hindered 
rotor part denoted by Hr(p, ip) and the analysis now fo- 
cuses on a one dimensional Hamiltonian due to the fact 
that P is a conserved quan tity. Note that this is con- 
sistent with the observation [108| by Lawton and Child 
regarding the evaluation of the frequency factor hui nm 
in eq. QT] The classical variables (P, p) are quantized as 
P = (ni + 712 + l)h = (m + l)h and p — (n\ — ri2)h = 
r?i;thereby the rotor barrier is different for different states 
|m, r) i.e., Vq — Vo(m,r). In this rotor representation, 
motion below and above the rotor barrier correspond to 
normal and local mode behaviour respectively. Exploit- 
ing the fact that Hr is a one dimensional Hamiltonian a 
semiclassical (WKB) expression for the local mode split- 
ting can be written down as 



A _ Wo -8 

2 _ 2ir 6 



(19) 



The frequency factor is loq = 2dHn/dp and the the tun- 
neling action integral 9 is taken between the two turn- 
ing points p(i/j*) = 0. Since the local modes are above 
the rotor barrier the turning points are purely imaginary 
ip* = ir/* = cosh _1 (i?^/Vb)/2 and thus 



dipp(tp) 



drj\J 8D(E R - Vq cosh(2r?)) (20) 



where, E R = fl 2 r 2 /8D and a crucial assumption has been 
made that Vo does not depend on p. Although such an 
asusmption is not strictly valid the estimates for the split- 
tings agreed fairly well with the exact quantum values. 

The analysis thus implicates a u>±:u>2 = 1:1 nonlinear 
resonance in the phase space for the dynamical tunneling 




t/T„ 



Figure 6: Three classes of energy flow dynamics, normal (N), 
local (L), and nonclassical (NC) illustrated for the coupled 
Morse system by computing the survival probabilities P(t) 
for the states |3, 2), [3, 0), and |3, 1} respectively. The cases 
L and NC are examples of dynamical tunneling whereas case 
N occurs classically. Note the difference in the timescales for 
complete energy transfer between modes for L and NC. In 
the case of NC significant probability builds up in the normal 
mode state |2, 2} (symbols). Very little probability buildup in 
the intermediate states |2, 1} (open circles) and |1,2) (filled 
circles) is seen in case L. 



between local modes and several interesting observations 
were made. First, the rotor Hamiltonian explains the 
coupling scheme outlined in Eq. (fTli)|) . Second, although 
there is a fairly strong coupling between |3, 0) and the 
intermediate state \2, 1) no substantial probability build- 
up was noticed in the state |2, 1) (cf. Fig. [6^)- On the 
other hand substantial probability does accumulate, as 
shown in Fig. [S] (NC), in the intermediate state |2,2) in- 
volved in the energy tran sfer process |3, 1) <-> |1,3). It 
was also commented [l04{ that t he d ouble well analogy 
provided by Davis and Heller [3, 115 1 was very different 
from the one that emerges from the hindered rotor analy- 
sis - local modes in the former case would be trapped be- 
low the barrier whereas they are above the barrier i n the 
latter case. Finally, a very important observation [lO^I 
was that small amounts of asymmetry between the two 
modes quenched the |3, 0) «-> |0, 3) process whereas there 
was little effect on the |3, 1) <-> |1,3) process. The dif- 
ferent effects of the asymmetry was one of the primary 
reasons for distinguishing between local and nonclassical 
states. 

Several questions arise at this stage. Is the rotor analy- 
sis applicable to the Lawton-Child and Davis-Heller sys- 
tems? Why is there a difference in the interpretation of 
the local modes, and hence the mechanism of dynami- 
cal tunneling, between the phase space and hindered ro- 
tor pictures? Is it reasonable to neglect the higher or- 
der uj\:u>2 — p:q resonances solely based on the relative 
strengths? If multiple resonances do exist then one has 
the possibility of their overlap leading to chaos. Is it still 
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possible to use the rotor formalism in the mixed phase 
space regimes? Some of the questions were answered in 
a work by Stefanski and Pollak wherein a periodic or- 
bit quantization technique was proposed to calcu late the 
splittings |l20| . Stefanski and Pollak pointed out [l20| an 
important difference between the Davis-Heller and the 
Sibert-Reinhardt-Hynes Hamiltonians in terms of the pe- 
riodic orbits at low energies. Nevertheless, they were able 
to show that a harmonic approximation to the tunneling 
action integral Eq. (1201) yields an expression for the split- 
ting which is identical to the one derived by assuming 
that the symmetric stretch periodic (unstable) orbit gives 
rise to a barrier separating the two local mode states - 
an interpretation that Davis and Heller provided in their 
workQ. Thus Stefanski and Pollak resolved an apparent 
paradox and emphasized the true phase space nature of 
dynamical tunneling. 

In any case it is worth noting that irrespective of the 
representation the key feature is the existence of a dy- 
namical barrier separating two qualitatively different mo- 
tions; the corresponding structure in the phase space had 
to do with the appearance of a resonance. Further sup- 
port for the importance of the resonance to tunneling 
between tori in the phase space came from a beautiful 
analysis by Ozorio de AlmeidaQ which, as seen later, 
forms an important basis for the recent advances. A dif- 
ferent viewpoint, using group theoretic arguments, based 
on the c once pt of dynamical symmetry breaking was 
advanced 121| by Kellman. Adapting Kellman's argu- 
ments, originally applied to the local mode spectrum of 
benzene, imagine placing one quantum in one of the lo- 
cal O-H stretching mode of H 2 0. Classially the energy 
remains localized in this bond and thus there is a low- 
ering or breaking of the C2 V point group symmetry of 
the molecule. However quantum mechanically, dynam- 
ical tunneling restores the broken symmetry. Invoking 
the more general permutation-inversion group |l22| and 
its feasible subgroups provides insights into the pattern 
of local mode splittings and hence insights into the en- 
ergy transfer pathways via dynamical tunneling. The 
dominance of a pathway of course cannot be established 
within a group theoretic approach alo ne a nd requires ad- 
ditional analysis. Note that recentlv |l23l | Babyuk, Wy- 
att, and Frederick have reexamined the dynamical tun- 
neling in the Davis-Heller and the coupled Morse sys- 
tems via the Bohmian approach to quantum mechanics. 
Analysing the relevant quantum trajectories Babyuk et 
al. discovered that there were several regions, at differ- 
ent times, wherein the potential energy exceeds the total 
energy (which includes the quantum potential). In this 
sense, locally, one has a picture that is similar to tunnel- 
ing through a one dimensional potential barrier. Inter- 
estingly such regions were associated with the so-called 
quasinodes which arise during the dynamical evolution 
of the density. Therefore the dynamical nature of the 
barriers is clearly evident, but more work is needed to 



understand the origin and distributions of the quasin- 
odes in a given system. Needless to say, correlating the 
nature of the quasinodes to the underlying phase space 
structures would be a useful endeavour. 



QUANTUM MECHANISM: VIBRATIONAL 
SUPEREXCHANGE 



The earlier works established the importance of the 
phase space perspective for dynamical tunneling using 
systems that possesed discrete symmetries and there- 
fore implicitly invoked an analogy to symmetric double 
well models. That is not to say that the earlier studies 
presumed that dynamical tunneling would only occur in 
symmetric systems. Indeed a careful study of the various 
papers reveal several insightful comments on asymmetric 
systems as well. However mechanistic details and quanti- 
tative estimates for dynamical tunneling rates were lack- 
ing. Important contributions in this regard were made 
by Stuchebrukhov , Meh ta, and Marcus nearly a decade 
ago [To! [To! EH, EH]. The experiments [aj that mo- 
tivated these studies have been discussed earlier in the 
introduction. In this section the key aspects of the mech- 
anism are highlighted. 

The inspiration comes from the well known superex- 
change mec hanism of long distance electron transfer 
in mol ecules Il27l Il28j . Stuchebrukhov and Marcus 
argued [103l Il05j that since the initial, localized (bright) 
state is not directly coupled by anharmonic resonances 
to the other zeroth-order states it is necessary to invoke 
off-resonant virtual states to explain the sluggish flow of 
energy. Specifically, it was noted that at any given time 
very little probability accumulates in the virtual states. 
In this sense the situation is very similar to the mech- 
anism proposed by Hutchinson, Sibert, and Hvnes 104 
as illustrated by Eq. (fT5)) . However Stuchebrukhov and 
Marcus exten ded t he mechanism, called as vibrational 
superexchange [l05j | . to explain the flow of energy between 
inequivalent bonds in the molecule and noted the surpris- 
ing accuracy despite the large number of virtual transi- 
tions invol ved in the process. As noted earlier, Hutchin- 
son's work |l06l | on state mixing in cyanoacetylene also il- 
lustrated the flow of energy between inequivalent modes. 
In order to illustrate the essential idea consider the hin- 
dered rotor Hamiltonian (cf. Eq. (fT8|) ): 



H R = ^ + V (P,p)cos(2iP) 



(21) 



where we have denoted / = AD /ft 2 . The free rotor ener- 
gies and the associated eigenstates are known to be 



E° R (k) 



21 



(22a) 
(22b) 
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The perturbation to the free rotor connects 
states that differ by two rotational quanta i.e., 
(k\V (P,p)cos(2i())\k') « V (P)/2 with \k - k'\ = 2. 
Clearly the local mode states \m,r), and \m, — r) with 
r 1 are not directly connected by the perturbation. 
Nevertheless the local mode states can be coupled 
through a sequence of intermediate virtual states with 
quantum numbers — (r — 2) < k < (r — 2). The effective 
coupling matrix element 



v (r_2) 

^ II 



Vo/2 



2 J-J- M(r) - E°Jl) 



(23) 



can be obtained via a standard application of high-order 
perturbation theory involving the sequence 



\m, r) — > \m, r — 2) 



|m,-(r-2)) - |? 



-r) 
(24) 

Note that the polyad quantum number m in the above se- 
quence is fixed and is consistent with the single resonance 
approximation. In addition Vo/2 -C £r(t) — En(r — 2) 
must be satisfied for the perturbation theory to be valid. 
Thus the splitting between the symmetry related local 
mode states |m,±r) is given by A m r = 2 T4ff . Interest- 
ingly, Stuchebrukhov and Marcus showed [l05j that V e e 
could be derived by analysing the semiclassical action 
integral 



dprfip) 



(25) 



with Tj = — iip and \p m \ = y/2I(E^(r) — Vq) being the 
minimum classical value of the momentum of the ro- 
tor excited above the barrier Vq. This suggests that the 
high-order perturbation theory, if valid, is the correct ap- 
proach to calculating dynamical tunneling splittings in 
multidimensions. An additional consequence is that the 
nonclassical mechanism of IVR is equivalent to dynam- 
ical tunneling as opposed to the earlier suggestion |l04j 
of an activated barrier crossing. As a cautionary note it 
must be mentioned that the above statements are based 
on insights afforded by near-integrable classical dynamics 
with two degrees of freedom. 

Although the example above corresponds to a sym- 
metric case the arguments are fairly general. In fact the 
vibrational superexchange mechanism is appropriate for 
describing quantum pathways for IVR in molecules. For 
example consider the acetylinic CH-stretch (vi) excita- 
tions in propyne, H3CCCH. Experiments by Lehmann, 
Scoles, and coworkers[92] indicate that the overtone state 
iv\ has a faster IVR rate as compared to the nearly 
degenerate v\ + 2vq combinatio n st ate. Similarly Go, 
Cronin, and Perry in their study [l24| found evidence for 
a larger number of perturbers for the 2v\ state than for 
the V\ + Vq state. The spectrum corresponding to both 
the first and the second overtone states implied a lack 



of direct low-order Fermi resonances. It was shown[126| 
that the slow IVR out of the 2v\ and 3^i states of propyne 
can be explained and understood via the vibrational su- 
perexchange mechanism. 

There is little doubt that the vibrational superex- 
change mechanism, as long as one is within the perturba- 
tive regimes, is applicable to fairly large molecules. How- 
ever several issues still remain unclear. The main issue 
has to do with the effective coupling between the initial 
bright state |b) = \v\,v\,--) and a nearly degenerate 
zeroth-order dark state |d) = \vf , > ■ • •) m the multidi- 
mensional state space. In the case of a lone anharmonic 
resonance there is only a single off-resonant sequence of 
states connecting |b) and |d) as in Eq. This se- 

quence or chain of states translates to a path in the state 
space. Note that the 'paths' being refered to here are 
nondynamical, although it might be possible to provide 
a dynamical interpretat ion b y analysing appropriate dis- 
crete action functional 103]. In general there are several 
anharmonic resonances of different orders and in such in- 
stances the number of state space paths can become very 
large and the issue of the relative importance of one path 
over the other arises. For instance it could very well be 
the case that |b) and |d) can be connected by a path 
with many steps involving low order resonances and also 
by a path involving few steps using higher order reso- 
nances. It is not a priori clear as to which path should 
be considered. An obvious choice is to use some pertur- 
bative criteria to decide between the many paths. Such 
a procedure was used by Stuch ebru khov and Marc us to 
create the tiers in their analvsis [l03j . More recently [l2i| 
Pearman and Gruebele have used the perturbative crite- 
ria to estimate the importance of direct high/low order 
couplings and low order coupling chains to the IVR dy- 
namics in the state space. Thus consider the following 
possible coupling chains, shown in Fig. [3J assuming that 
the states |b) and |d) are separated by n quanta in the 
state space i.e., \ v t ~ v t \ = 71 



|ix> ^ |d) 

lix) ^ |i 2 ) |d) 



Note that n is the distance between |b) and |d) in the 
state space and thus identical to Q in Eq. (fTUjl . In the 
above the first equation indicates a direct coupling be- 
tween the states of interest by an anharmonic resonance 

(n) 

coupling Vfr d of order n. In the second case the coupling 
is mediated by one intermediate state |ii) involving two 
anharmonic resonances V^™ 1 ^ and V^f 1 with n = 774+77,2. 
Each step of a given state space path, coupling two states 
|i) and |j) at order m, is weighted by the perturbative 
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term 



C 



(m) 



-1/2 



(m) 



(26) 



(m) 



for A4»^ m) 

(m) 



1 



for A^.«^ 



with AM 



Ef - E° and V>™ } = (i|y (m) |j). A specific 
path in the state space is then associated with the prod- 
uct of the weightings for each step along the path. For ex- 
ample the state space path above involving two interme- 
diate states is associated with the term C^^C^^C^^ . 
It is not hard to see that such products of C correspond to 
the terms contributing to the effective coupling as shown 
in Eq. (fT6)) and Eq. (|23|) . The various terms also con- 
tribute to the effective coupling ipQ (cf. Eq. (flOf ) in the 
Leitner-Wolynes criteria[87| for the transition from local- 
ized to extended states. 

The crucial issue is wether or not one can identify 
dominant paths, equivalently the key intermediate states, 
based on the perturbative criteria alone. In a rigorous 
sense the answer is negative since such a criteria ignores 
the dynamics. Complications can also arise due to the 
fact tha t eac h segment of the path contributes with a 
± phase [123 • Furthermore one would like to construct 
the explicit dynamical barriers in terms of the molecu- 
lar parameters and conserved or quasi-conserved quan- 
tities. Th e obs ervation that the superexchange V c r can 
be derived 103] from a semiclassical action integral pro- 
vides a clue to some of the issues. However an explicit 
demonstration of such a correspondence has been pro- 
vided only in the single resonance (hindered rotor) case. 
In the multidimensional case the multiplicity of paths ob- 
scures this correspondence. The next few sections high- 
light the recent advances which show that clear answers 
to the various questions raised above come from view- 
ing the phenomenon of dynamical tunneling in the most 
natural representation - the classical phase space. 



words 'remote' and 'distant' are not immediately clear 
and fraught with conceptual difficulties in a multidimen- 
sional phase space setting. In essence Heller's statement, 
more appropriately called as a conjecture, is an effort to 
provide a phase space picture of the superexchange mix- 
ing between two or more widely separated zeroth-order 
states in the state space. Some of the examples in the 
later part of the present review, hopefully, demonstrate 
that the conjecture is reasonable. Interestingly though in 
the same review it is mentioned fl3| that in the presence 
of widespread chaos the issue of tunneling as a mecha- 
nism for IVR is of doubtful utility. Similar sentiments 
were echoed by Stuchebrukhov et al. who, in the case of 
high dimensional phase spaces, envisaged partial ruptur- 
ing of the invariant tori and hence domination of IVR by 
chaotic diffusion rather than dynamical tunneling. This 
too has to be considered as a conjecture at the present 
time since the timescales and the competition between 
chaotic diffusion (classically allowed) and dynamical tun- 
neling (classically forbidden) in systems with three or 
more degrees of freedom has not been studied in any de- 
tail. The issues involved are subtle and extrapolating the 
insights from studies on two degrees of freedom to higher 
degrees of freedom is incorrect. For instance, in three and 
higher degrees of freedom the invariant tori, ruptured or 
not, do not have the right dimensionality to act as bar- 
riers to inhibit diffusion. Therefore, is it possible that 
classical chaos could assi s t or i nhib i t dynami c al tunnel- 
ing? Several studies [TH, EM EH, \±M, EU EH EH 



1 168L Il69l | carried out during the early nineties 
till the present time point to an important role played 
by chaos in the process of dynamical tunneling. We start 
with a brief review of these studies followed by very de- 
tailed investigations into the role played by the nonlinear 
resonances. Finally, some of the very recent work is high- 
lighted which suggest that the combination of resonances 
and chaos in the phase space can yet play an important 
role in IVR. 



Chaos-assisted tunneling (CAT) 



DYNAMICAL TUNNELING: RECENT WORK 

Nearly a decade ago Heller wrote a brief review [l]| 
on dynamical tunneling and its spectral consequences 
which was motivated in la rge part by the work of Stuche- 
brukhov and Marcus |l03l |. Focusing on systems with two 
degrees of freedom in the near-integrable regimes, charac- 
teristic of polyatomic molecules at low energies, Heller ar- 
gued that the lO^-lO -2 cm broadening of lines is due 
to dynamical tunneling between remote regions of phase 
space facilitated by distant resonances. This is an in- 
triguing statement which could perhaps be interpreted 
in many different ways. Moreover the meaning of the 



One of the first studies on the possible influen ce of 
chaos on tunneling was made by Lin and Ballentine [l30l 



13U |133| | . These authors investigated the tunneling dy- 
namics of a driven double well system described by the 
Hamiltonian 

H (x, p, t) = -^—p 2 + bx 4 - dx 2 + X lX cos(cjt) (27) 
2m 

which has the discrete symmetry H(x,p,t) = 
H(~x, —p, t + To/2) with To = 2tt/u>. In the absence of 
the monochromatic field (Ai = 0) the system is integrable 
and the phase space is identical to that of a one dimen- 
sional symmetric double well potential. However with 
Ai 7^ the system is nonintegrable and the phase space 
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is mixed regular-chaotic. Despite the mixed phase space 
the discrete symmetry implies that any regular islands in 
the phase space will occur in symmetry related pairs and 
the quantum floquet states will occur as doublets with 
even/odd symmetries. A coherent state (wavepacket) lo- 
calized on one of the regular islands will tunnel to the 
other, clasically disconnected, symmetry related island. 
Thus this is an example of dynamical tunneling. An im- 
portant observation by Lin and Ballentine was that the 
tunneling was enhanced [l3Ct[l3l| by orders of magnitude 
in regimes wherein significant chaos existed in the phase 
space. For example with Ai = the ground state tun- 
neling time is about 10 6 To whereas with Ai = 10 the 
tunneling time is only about 10 2 Xo- Futhermore strong 
fluctuations in the tunneling times were observed. In 
Fig. [7] a typical example of the fluctuations over several 
orders of magnitude are shown. The crucial thing to 
note, however, is that the gross phase space structures 
are similar over the entire range shown in Fig. [7l Thus 
although the chaos seems to influence the tunneling the 
mechanistic insights were lacking i.e., the precise role of 
the stochasticity was not understood. 




Figure 7: Tunneling time for the Lin-Ballentine system as 
a function of the field strength Ai. The parameters are 
chosen 130] as m — 1, b = 0.5, and d — 10. The phase 
space for Ai = 7.8 in the upper left panel shows the two reg- 
ular islands embedded in the chaotic sea. The initial state 
corresponds to a coherent state localized on the left regular 
island (indicated by an arrow). The three phase space strobe 
plots at increasing values of Ai indicated in (a) show the local 
structure near the left island. Despite the phase space being 
very similar the tunneling times fluctuate over many orders 
of magnitude. Figure courtesy of Astha Sethi (unpublished 
work) . 



mann et al. who showed that it is possible to suppress 
the tunneling in the driven double well by an appro- 
priate choice of the field parameters (Ai,w). The ex- 
planation for such a coherent destruction of tunneling 
and the fluctuations comes from analysing the the flo- 
quet lev el mo tions with Ai. Gomez Llorcnte and Plata 
provided 135] perturbative estimates within a two-level 
approximation for the enhancement/suppression of tun- 
neling. A little later Utermann, Dittrich, and Hanggi 
showed [l36| | that there is a strong correlation between 
the splittings of the floquet states and their overlaps 
with the chao tic part of the phase space. Breuer and 
Holthaus 214j highlighted the role of classical phase space 
structures in the d riven double well system . Subse- 
quently other works [m EM EM, EM EM OH EH 
on a wide variety of driven systems established the sensi- 
tivity of tunneling to the classical stochasticity. An early 
discussion of dynamical tunneling a nd t he influence of 
dissipation can be found in the work|l44j] by Grobe and 
Haake on kicked tops. A comprehensive account of the 
various studies on driven systems can be found in the 
review [l45| by Grifoni and Hanggi. Such studies have 
provided, in recent times, important insights into the 
proces s of c oherent control of quantum processes. Recent 
review [l46l | by Gong and Brumer, for instance, discusses 
the issue of quantum control of classically chaotic sys- 
tems in detail. Perhaps it is apt to highlight a historical 
fact mentioned by Gong and Brumer - coherent control 
emerged from two research groups engaged in studies of 
chaotic dynamics. 

In the Lin-Ballentine example above the perturbation 
(applied field) not only increases the size of the chaotic 
region but also affects the dynamics of the unperturbed 
tunneling doublet. Hence the enhanced tunneling, rel- 
ative to the unperturbed or the weakly perturbed case, 
cannot be immediately ascribed to the increased amount 
of chaos in the phase space. To observe a more direct 
influence of chaos on tunneling it is necessary that the 
classical phase space dynamics scales with energy. Inves- 
tigations of the coupled qua rtic oscilla tor system by Bo- 
higas, Tomsovic, and Ullmo [147l Il48j | provided the first 
detailed view of the CAT process. The choice of the 
Hamiltonian 



H(p,q) = ±p 2 + a(X) 



|+6, 2 4 + 2A^l 



(28) 



Important insights came from the work 134j by Gross- 



with b — 7r/4 was made in order to study the classical- 
quantum correspondence in detail. This is due to the 
fact that the potential is homogeneous and hence the dy- 
namics at a specific energy E is related to the dynamics 
at a different energy by a simple rescaling. Consequently 
one can fix the energy and investigate the effect of the 
semiclassical limit fi — >■ on the tunneling process. The 
classical dynamics is integrable for A = and strongly 
chaotic for A = — 1. Again, note that the Hamiltonian 
has the discrete symmetry q — > q, p — > — p and thus the 
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quantum eigenstates are expected to show up as symme- 
try doublets. In contrast to the driven system discussed 
above, and similar to the Davis-Heller system discussed 
earlier, the system in eq. ([25]) does not have any barriers 
arising from the two dimensional quartic potential. The 
doublets therefore must correspond to quantum states 
localized in the symmetry related islands in phase space 
shown in Fig. [8] Thus denoting the states localized in 
the top and bottom regular islands by \T) and \B) re- 
spectively the doublets correspond to the usual linear 
combinations |±) = (|T) ± \B))/V2. 




Figure 8: In (a) the q2 — surface of section is shown for 
E = 1 and the parameter values (A, b) = (— 0.25,7r/4) for 
the Hamiltonian in Eq. (|28p . The large islands are part of a 
1:1 resonance island chain and are related to each other by 
symmetry, (b) Details of the upper island and (c) shows the 
configuration space plot of the symmetrically related tori. 



Bohigas, Tomsovic, and Ullmo focused on quantum 
states associated with the large regular islands seen in 
Fig. [S] Specifically, they studied the variations in the 
splitting AE = \E + — E-\ associated with the doublets 
|±) upon changing fi and the coupling A. Since the states 
are residing in a regular region of the phase space it is 
expected, based on one dimensional tunneling theories, 
that the splittings will scale exponentially with H i.e 



AE ~ exp(— 1/h). However the splittings exhibit [147 1 
fluctuations of several orders of magnitude and no obvi- 
ous indication of a predictable dependence on h seemed 
to exist. Note that a similar feature is seen in Fig.[7]show- 
ing the tunneling time fluctuations of the drive n double 
well system. More precisely, it was observed |l47l | that for 
the range of A that corresponded to integrable or near- 
integrable phase space AE roughly follows the exponen- 
tial scaling. On the other hand for A = A* w —0.2 chaos 
becomes pervasive and AE exhi bits severe fluctuations. 
Tomsovic and Ullmo argued[148] that the fluctuations in 
AE could be traced to the crossing of the quasidegen- 



erate doublets with a third irregular state. By irregu- 
lar it is meant that when viewed from the phase space 
the state density, represented by either the Wigner or 
the Husimi distribution function, is appreciable in the 
chaotic sea. The irregular state, in contrast to the regular 
states, do not come as doublets since there is no dynam- 
ical partition of the chaotic region into mutually exclu- 
sive symmetric parts. Nevertheless the irregular states do 
have fixed parities with respect to the reflection opera- 
tions. Thus chaos-ass isted tunn eling is necessarily atleast 
a three-state process [l48[[l4^ . In other words, assum- 
ing that the even-parity irregular state |C) with energy 
E c couples to |+) with strength v a relevant, but min- 
imal, model Hamiltonian in the symmetrized basis has 
the form: 



E + e 
If = | E-e v 

v E r 



(29) 



with e being a direct coupling between \T) and \B). If e 
is dominant as compared to v then one has the usual two 
level scenario. However if « > e rj then the splitting 
can be approximated as 



AE 



\v\ E-E c <^v. 



(30) 



Hence varying a parameter, for example A, one expects 
a peak of height \v \ as E c cr osses the E+ level. Detailed 
theoretical studies [15Cj, Il51j of the dynamical tunneling 



process in annular billiards by Frischat and Doron con- 
firmed the crucial role of the classical phase space struc- 
ture. Specifically it was found that the tunneling be- 
tween two symmetry related whisphering gallery modes 
(corresponding to clockwise and counterclockwise rotat- 
ing motion in the billiard) is strongly influenced by quan- 
tum states that inhabit the regular-chaotic border in the 
phase space. Such states were termed as "beach states" 
by Frischat and Doron. Soon thereafter Dembowski et 
al. provided experimental evidence [l6| for CAT in a mi- 
crowave annular billiard simulated by means of a two- 
dimensional electromagnetic microwave resonator. De- 
tails of the experiment can be found in a recent paper [l7l] 
by Hofferbert et al. wherein support for the crucial role 
of the beach region is provided. It is, however, signif- 
icant to note that the regular-chaotic border in case of 
the annular billiards is quite sharp. 

An important point to note is that despite the intu- 
tively appealing three-state model, determining the cou- 
pling v between the regular and chaotic states is non- 
trivial. This, in turn, is related to the fact that accu- 
rate determination of the positions of the chaotic levels 
can be a difficult task and determining the nature of a 
chaotic state for systems with mixed phase spaces is still 
an open p roble m. Such difficulties prompted Tomsovic 
and Ullmo[148[ to adopt a statistical approach, based 
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on random matrix theory, to determine the tunneling 
splitting distribution P(AE) in terms of the variance 
of the regular-chaotic coupl ing i.e., v 2 . In a later work 
Leyvraz and Ullmo showed [152j that P(AE) is a trun- 
cated Cauchy distribution 



P(AE) = I- XE 



7T (AE) 2 + (AE) 2 



(31) 



with AE being the mea n splitting. On the other hand 
Creagh and Whelan 154l 1 1 561 ] showed that the splitting 
between states localized in the chaotic regions of the 
phase space but separated by an energetic barrier is char- 
acterized by a specific tunneling orbit. Specifically, the 
resulting splitting distribution depends on the stability 
of th e tun neling orbit and is not universal. In certain 
limits 154j the splittings obey the Porter-Thomas distri- 
bution. Thus the fluctuations in AE in the CAT pro- 
cess are different from those found in the usual double- 
well system. Note that in case of CAT the distribu- 
tion pertains to the splitting between states localized 
in the regular regions of phase space embedded in the 
chaotic sea. The reader is r eferre d to the excellent dis- 
cussion by Creagh for details []J57j]. Nevert heles s, in a re- 
cent work Mouchet and Delande observed [l58| that the 
P(AE) exhibits a truncated Cauchy behaviour despite 
a near-integrable phase space. Similar observations have 
been made by Carvalho and Mijolaro in their study of 
the splitting distribution in the annular billiard system 
as well |l59j | . This indicates that the Leyvraz-Ullmo dis- 
tribution Eq. (|3"Tj) is not sufficient to characterize CAT. 
A first step towards obtaining a semiclassical estimate of 
the regular -chao tic coupling was taken by Podolskiy and 
Narimanov |l60l | . Assuming a regular island seperated by 
a structureless, on the scale of ft, chaotic sea they were 
able to show that: 



AE cx ft 



T(A, 2A) 
T(A + 1,0) 



(32) 



with A = A/nft and a system specific proportionality 
factor which is independent of ft. The phase space area 
covered by the regular island is denoted by A. Applica- 
tion of the theory by Podolskiy and Narimanov to the 
splittings between near-degenerate optical modes, local- 
ized on a pair of symmetric regular islands in phase space, 
in a non-integrable microc avity yielded very good agree- 
ment with the exact data[l60(. In fact Podolskiy and 
Narimanov have recently shown 16l| that the lifetimes 
and emission patterns of optical modes in asymmetric 
microresonators are strongly affected by CAT. Encour- 
aging results were obtained in the case of tunneling in 
periodically modulated optial lattices as well. However 
the agreement in the splittings displayed significant de- 
viations in the deep semiclassical regime i.e., large vaues 
of ftT . Interestingly the critical value of ft -1 beyond 
which the disagreement occurs correlates with the exis- 
tence of plateau regions, discussed in the next section (cf. 



Fig-HU), m the plot of log(A_E) versu s ft -1 . Such plateau 
regions have been noted earlier 1 1621 by Rqncag lia et al. 
and in several other recent studiesll78l Il82l \l83k as well. 



The qualitative picture that emerged from the numer- 
ous studies is that CAT process is a result of competition 
between classically allowed (mixing or transport in the 
chaotic sea) and classically forbidden (tunneling or cou- 
pling between regular regions and the chaotic sea) dy- 
namics. Despite the significant advancements, the mech- 
anism by which a state localized in the regular island 
couples to the chaotic sea continued to puzzle the re- 
searchers. It might therefore come as a surprise that 
based on recent studies there is growing evidence that 
the explanation for the regular-chaotic coupling lies in 
the nonlinear resonances and partial barriers in the phase 
space. The work by Podolskiy and Narimanov provides 
one possible answer but an equally strong clue is hidden 
in the plateaus that occur in a typical log(A-E) versus ft -1 
plot. It is perhaps reasonable to expect that the theory 
of Podolskiy and Narimanov is correct when phase space 
structures like the nonlinear resonances and partial bar- 
riers, ignored in the derivation of Eq. (|32p . in the vicin- 
ity of regular-chaotic border regions are much smaller in 
area as compared to the Planck constant ft. However, 
for sufficiently small ft the rich hierarchy of structures 
like nonlinear resonances, partially broken tori, cantori, 
and even vague tori in the regular-chaotic border region 
have to be taken into consideration. It is significant to 
note that Tomsovic and Ullmo had already recognized 
the importance of accounting for mechanisms which tend 
to limit classical transport [1481 ] and proposed generalized 
random matrix ensembles as a poss i ble a pproach. There 
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166, 16 



also exist a number of studies 163, 
involving less than three degrees of freedom, which high- 
light the nature of the regular-chaotic border and their 
importance to tunneling between KAM tori. 

A different perspective on the influence of chaos on dy- 
namical tunneling came from the detailed investigations 
by Ikeda and coworkers involving the studv |l68 . 169l 17Ct 



ml Hz! [nil of semiclassical propagators in the complex 



phase space. In particular, deep insights into CAT were 
gained by examining the so called "Laputa chains" which 
contribu te domin antly to the propagator in the pre sence 
of chaos [l68i Il69j . For a detailed review the paper [l70j ] 
by Shudo and I keda is highly reccommended. Hashimoto 
and TnkatMika 17."> discuss a situation wherein dynam- 
ical tunneling leads to mixing between states localized 
about unstable fixed points in the phase space. However 
in this case the transport is energetically as well as dy- 
namically allowed and it is not very clear whether the 
term 'tunneling' is an appropriate choice. 

In the following subsections recent progress towards 
quantitative prediction of the average tunneling rates, 
and hence the coupling strengths between the regular 
and chaotic states, are described. The key ingredient to 
the success of the theory are the nonlinear resonances in 
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the phase space. In fact, depending on the relative size of 
h, one needs to take multiple nonlinear resonances into 
account for a correct description of dynamical tunneling. 
If one asociates the h variation with the related density 
of states in a molecular system then it is tempting to 
claim that such a mechanism must be the correct semi- 
classical limit of the Stuchebrukhov-Marcus vibrational 
superexchange theory. The discussions in the next sec- 
tion indicate that such a claim is indeed reasonable. 



Resonance-assisted tunneling (RAT) 

The early investigations of dynamical tunneling made 
it clear that nonlinear resonances played an important 
role in the phenomenon of dynamical tunneling. How- 
ever in the majority of the studies one had a single reso- 
nance of low order dominating the tunneling process and 
hence an effectively integrable classical dynamics. Very 
few attempts were made for systems which had multiple 
resonances and thus capable of exhibiting near-integrable 
to mixed phase space dynamics. Moreover studying the 
tunneling process with varying Ti is necessary to implicate 
the nonlinear resonances and other phase space struc- 
tures without ambiguity. Definitive answers in this re- 
gard were provi ded by Brodier, Schlagheck, and Ullmo in 
their study [l76j of time-dependent Hamiltonians. In this 
section we provide a brief outline of their work and refer 
the reader to some of the recent reviews 177, llTSL Il79l | 
which provide details of the theory and applications to 
several other systems. 

In order to highlight the salient features of RAT con- 
sider a one degree of freedom system that evolves un- 
der a periodic time-dependent Hamiltonian H(q,p,t) = 
H(q,p,t + t). The phase space structure in this case 
can be easily visualized by constructing the stroboscopic 
Poincare section i.e., plotting (q,p) at integer values of 
the period t = nr. We are interested in the generic 
case wherein the phase space exhibits two symmetry- 
related regular islands that are embedded in the phase 
space. Typically the phase space has various nonlinear 
resonances that arise due to frequency commensurability 
between the external driving field frequency lu = 2it/t 
and the system frequencies. Fig. [9] shows an example of 
a near-integrable phase space wherein a few resonances 
are visible. For the moment assume that a prominent 
r:s, say the 10:1 resonance in Fig.[9j nonlinear resonance 
manifests in the phase space. The motion in the vicin- 
ity of this r:s resonance can be analysed using secular 
perturbation theory. In order to do this one introduces 
a time-independent Hamiltonian Ho(q,p) that approx- 
imately describes the regular motion in the islands. In 
particular assume that appropriate action-angle variables 
(1,6) can be introduced such that H (q,p) — > Hq(I) 
and thus the total Hamiltonian can be expressed as 
H(I,6,t) = H (I) + V(I,6,t). The term V represents 



a weak perturbation in the center of the island. A r:s 
nonlinear resonance occurs when 



with f2 r ., = 



dH Q 



dl 



(33) 



The reso nant action is denoted by I r:s . Using 
techniques[180] from secular perturbation theory one can 
show that the dynamics in the vicinity of the nonlinear 
resonance is described by the new Hamiltonian [l80l | 



(I 
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V(i,#,t) 



(34) 



with V(J,0,t) = V(I,$ + Q r:s t,t), and the "effective 
mass" parameter m r:s = (d 2 Ho/dI 2 )^ 1 (I r:s ) is related 
to the anharmonicity of the system. The new angle 
d = 6 — Q r:s t varies slowly in time near the resonance 
and hence, according to adiabatic perturbation theory, 
one can replace V by its time average over r-periods of 
the driving field 



V{I,0) = — [ dtV(L 
rT Jo 



0,t) 



(35) 



Furthermore, writing the perturbation in terms of a 
Fourier series and neglecting the action dependence of 
the Fourier coefficients i.e., evaluated at I = I r:s the 
effective integrable Hamiltonian is obtained as: 



(I - I ) 2 °° 
Hr:.(I,0) = 9 +^Ffc(I r: s)cos(fcn9+0 fc ) (36) 

k=l 



The above effective Hamiltonian couples the zeroth or- 
der states |n) and \n + kr) with a matrix element 
14 exp(i(/)fc)/2. Here the quantum numbers refer to the 
zeroth-order KAM tori inside the large regular island. 
Thus, the eigenstates of H r:s are admixtures of unper- 
turbed states obeying the "selection rule" \n' — n\ = kr 
and perturbatively given by 



= In) 



y Vk 
^ AEn+kr 



n 



- kr) 



V — 

kW^O ^En+krAE n+ kr+k'r 



\n + kr + k'r) (37) 



where = Vfcexp(i^fc) and AEj = E n — Ej. Associ- 
ating the quantized actions = Ti(n + 1/2) with states 
localized in the regular islands the energy difference 



E n — E n ' = 



1 



-(In ~ In>)(In + In> ~ 2Ir:s) (38) 



2m r 

becomes small if /„ + I n ' ~ 2/ r:s . In other words, the 
states \n) and \n') are coupled strongly if they straddle 
the r:s resonance symmetrically. Consequently the non- 
linear resonance provides an efficient route to couple the 
ground state with Iq < I r:s to a highly excited state with 

Ikr ^ Ir:s- 
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Once again a crucial issue has to do with the the dom- 
inant pathways that couple |n) and \n + kr). Clearly the 
two states can couple in a one-step process via or with 
a fc-step process via V%. The analysis by Brodier et at, 
realizing the rapid decrease of the Fourier terms Vk with 
k, reveals |l77| that the /c-step process is dominant in the 
semiclassical limit 7i — > 0. Note that within the effec- 
tive Hamiltonian H r[s one assumes the dominance of a 
single nonlinear resonance and hence the focus is on the 
r:s resonance and its higher harmonics. Neglecting the 
Fourier components with k > 1 and denoting 2V r:s = V\ 
one obtains the standard pendulum Hamiltonian 



{I - I r 



2m r 



+ 2V r:s cos(n?) (39) 



For periodically driven systems the floquet states can be 
obtained by diagonalizing the one-period time evolution 
operator. The near degenerate floquet states are charac- 
terized by the splitting 



A0„ 



— u+ 



(40) 



between the quasienergies, or eigenphases, of the sym- 
metric and antisymmetric state. Within the perturba- 
tion scheme the ground state (n — 0) splitting is then 
obtained as 



Ha 



(r:s), 
kr 



] A(j) 



(0) 
At 



with 



A 



(r:s) — 
kr 



\k\ 

(kr\0) = n 



V r:s 



i=i 



En — E, 



Ir 



(41) 



(42) 



being a measure of the contribution of the fcr-th ex- 
cited state to the perturbed ground state. The quantities 
are the unperturbed splittings resulting from the 
integrable limit Hamiltonia n for which extensive semi- 
classical insights exist 0, [ll], Il8lj . Note that the above 
expression for is essentially the same as the one 

obtained by Stuchebrukhov and Marc us within their vi- 
brational superexchange approach 103| . Here one has the 
generalization for driven systems whose phase space ex- 
hibits a general r:s nonlin ear resona nce. 

Brodier et al. illustrated 176l 177 1 the RAT mechanism 
beautifully by studying the kicked Harper model 

oo 

H(p, q, t) = cos(p) + r ^2 nT ) cos(g) (43) 

n— — oo 

whose dynamics is equivalent to the symplectic map 

Pn+l = Pn + Tsm(q n ) 

q n +i = q n -Tsin(p„+i) (44) 

which describes the evolution of (p, q) from time t = nr 
to t = (n + 1)t. For r = 1 the phase space, shown in 
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Figure 9: The classical phase space (upper panel) and quan- 
tum eigenphase splittings (lower panel) of the kicked Harper 
system with r = 1. The splittings are calculated for the n-th 
excited state at TV = 2ir/h = 6(2n + 1). Variation with JV is 
equivalent to fixing the state localized on the torus with action 
I — tv / 6 and varying %. Exact quantum splittings (circles) are 
well reproduced by semiclassical prediction (thick line) based 
on the 16:2, the 10:1, and the 14:1 resonances. The dashed 
line is the integrable estimate obtained using the Hamiltonian 
in Eq. (|45[) and clearly inaccurate. Note that multiple pre- 
cision arithmetics was used to compute eigenphase splittings 
below the ordinary machine precision limit. Figure courtesy 
of P. Schlagheck. 



Fig. [HI is near-integrable and exhibits a large regular is- 
land around (p, q) — (0,0). The 2n periodicity of the sys- 
tem in both p and q implies the existence of a symmetry 
related island at (p, q) = (0, 2ir). The floquet states thus 
come in nearly degenerate pairs with eigenphase split- 
tings A(j) n associated with the n-th excited state in the 
regular island. Since the system is near-integrable for 
t = 1 it is possible to use classical perturbation theory 
to construct a time-independent Hamiltonian 



Ha(p,q) = cosp + cos q — — smp sing 



— (cos p sin 2 q + cos q sin 2 p) 
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sin 2p sin 2q 



(45) 



which is a very good integrable approximation to the 
kicked Harper map. Indeed the phase space for the exact 
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time-dependent Hamiltonian for t — 1 and the integrable 
approximation in Eq. (|45[) are nearly identical |l76| . The 
crucial difference, however, is that the integrable approx- 
imation does not have the various nonlinear resonances 
which are visible in Fig. [HI The consequence of negelct- 
ing the nonlinear resonances are catastrophic as far as 
the splittings are concerned and from Fig. [9] it is clear 
that any estimate based on the integrable approximation 
is doomed to fail. 

Such pendulum-like Hamiltonians have been used in 
the very early studies of dynamical tunneling and were 
also central to the ideas and conjectures formulated by 
Hellerpjl in the context of IVR and spectral signatures. 
However the importance of the contribution by Brodier, 
Schlagheck, and Ullmo lies in the fact that a rigorous 
semiclassical basis for the RAT mechanism was estab- 
lished. Although the detailed study was in the context 
of a one-dimensional driven system, the insights into the 
mechanism are expected to be valid in larger class of 
systems. These include systems with two degrees of free- 
dom and perhaps to some extent systems with three or 
more degrees of freedom as well (see discussions in the 
introduction and the last section). For instance it can 
be shown that the RAT mechanism will always dominate 
the regular tunneling process in the semiclassical limit. 
Furthermore, discussions on the importance of the higher 
harmonics of a given nonlinear resonance and the criteria 
for identifying and including multiple resonances in the 
path between \n) and \n') are established. In this regard 
an important point made by Brodier et al. is that the 
condition for including a r:s resonance into the coupling 
path is not directly related to the size of the resonance is- 
lands - a crucial point that emerges only upon investigat- 
ing the analytic cont inua tion of the classical dynamics to 
complex phase spacefill}. Moreover, provided the obser- 
vations are general enough, the arguments could possibly 
shed some light on the issue of the competition between 
low and high orde r an harmonic resonances for IVR in 
large molecules 87 , 129l | . 



RAT precedes CAT? 

Up until now the two processes of CAT and RAT 
have been discussed separately. Numerous model studies 
showed that the nonlinear resonances are important for 
describing the dynamical tunneling process in the near- 
integrable regimes. At the same time specific statisti- 
cal signatures could be associated with the influence of 
chaos on dynamical tunneling. One sticky issue still re- 
mained - the mechanism of the coupling between a reg- 
ular state and the chaotic sea. As noted earlier, the 
assumption of a structurless chaos-regular border and 
hence their unimportance for CAT process typi cally gets 
violated for sufficiently small h. Many authors [147], [l51| 
suspected that the nonlinear resonances could play an 
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Figure 10: Sketch of the effective Hamiltonian matrix that 
describes tunneling between the symmetric quasi-modes in 
the two separate regular islands. The regular parts (upper 
left and lower right band) includes only components that are 
coupled to the island's ground state by the r:s resonance. 
The chaotic part (central square) consists of a full sub-block 
with equally strong couplings between all basis states with 
actions beyond the outermost inva riant torus of the islands. 
Reproduced horn Schlagheck et al. [l78l |. 



important role in the mixed phase space scenario as well. 
However the precise role of the nonlinear resonances in 
determining the regular-chaotic coupling was not clear. 
An important hrst step in this dire ction was recently 
taken by Eltschka and Schlagheck [l82|. By extending the 
resonance-assisted mechanism Eltschka and Schlagheck 
argue that the average splittings i.e., ignoring the fluc- 
tuations associated with multistate avoided crosings, can 
be determined without any adjustable parameters. 

In order to illustrate the idea note that in the mixed 
regular-chaotic regime invariant tori corresponding to the 
regular region are embedded into the chaotic sea. A typ- 
ical case is shown in Fig. [IT] for the kicked Harper map 
with r = 2. Therefore only a finite number of invariant 
tori are supported by the regular island with the outer- 
most torus corresponding to the action I c . Assuming the 
existence of a dominant r:s nonlinear resonance implies 
that the ground state in the regular region |0) is coupled 
efficiently to the unperturbed state \kr). If \kr) is lo- 
cated beyond the outermost invariant torus of the island 
i.e., I(k—i)r < Ic < Ikr then one expects a significant 
coupling to the unperturbed states located in the chaotic 
sea. Within this scheme the effective Hamiltonian can be 
thought of as having the form shown in Fig. [TOl 

The effective coupling between the ground state and 
the chaos obtained by eliminating the intermediate states 
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within the regular island is given by 

k— 1 

V eS = V r:s I] - Vr:s (46) 

and follows from the near-integrable resonance-assisted 
theory. The chaotic part of the effective Hamiltonian 
is modeled[l82| by a random hermitian matrix from the 
Gaussian orthogonal ensemble. In particular, assuming 
the validity of the truncated Cauchy distribution for the 
splittings, the mean eigenphase splittings for periodically 
driven systems assumes the simple form 

<A^oH(^) 2 (47) 

with a universal prediction for the logarithmic variance 
of the splittings 

([ln(A0o)-(ln(A0 o ))] 2 )-^ (48) 

In other words the actual splittings might be enhanced or 
suppressed as compared to the mean by about exp(7r/2) 
independent of the value of fi and other external param- 
eters. In Fig. [ITJ the expression for the mean splitting 
in Eq. (|47f is compared to the exact quantum results in 
the case of the kicked Harper at r = 2. One immedi- 
ately observes that the agreement is relatively good and 
the approximate plateau like structures are also repro- 
duced. The plateaus are related to the number of per- 
turbative steps in Eq. (|46[) that are required to connect 
the ground state with the chaotic domain. In order to 
understand the plateaus note that increasing N in Fig. [Til 
corresponds to decreasing fi since the state is fixed. As fi 
decreases there comes a point when the state \kr) local- 
izes on, rather than being located beyond, the outermost 
invariant torus of the regular island. Thus the coupling 
to the chaotic domain via Eq. ([46|) now requires (k + 1) 
steps rather than k steps. 

The mechanism for CAT just outlined has been con- 
firmed for a variety of systems. Eltschka and Schlaghcck 
confirmed that the plateau structures appear for the 
kicked rotor c ase as well. Wimberger et al. very re- 
cently showed |l83l | that the decay of the nondispersive 
electronic wavepackets in driven ionizing Rydberg sys- 
tems is also controlled by the RAT mechanism. Mouchet 
et al. performed a fairly detailed analysis of an effective 
Hamiltonian that can be realized in experiments on cold 
atoms and provide clear evidence for the important role 
played by th e var ious nonlinear resonances [l84j. A ex- 
tensive study [l85l ] . by Sheinman et al., of the lifetimes 
associated with the decay of the quantum accelerator 
modes revealed regimes wherein the RAT mechanism is 
crucial. In this context it is interesting to ask wether the 
mechanism would also explain the CAT process in the 
Bohigas-Tomsovic-Ullmo system described by the Hamil- 
tonian in Eq. (j2"8"|) . It is evident from the inset in Fig. [SJa 
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Figure 11: (Colour online) A typical mixed regular-chaotic 
phase space obtained for the kicked Harper map with r = 2 
is shown in the upper panel. A prominent 8:2 nonlinear res- 
onance is visible. The lower panel shows the quantum eigen- 
phase splittings calculated for the ground state within the 
regular island as a function of TV = 2n/h. The thick solid 
line (red) is the semiclassical prediction based on the 8:2 res- 
onance for the average splittings along with the logarithimic 
standard deviations as dashed lines (red) . Th e long dashed 
curve (green) is the Podolskiv-Narimanov [l60| prediction for 
the same system. Figure courtesy of P. Schlagheck. 



that nonlinear resonances do show up around the chaos- 
regular border but a systematic study of the original sys- 
tem from the recent viewpoint s is y et to be done. Similar 
questions are being addressed [186J in the Lin-Ballentine 
system (cf. Eq. ([57])) i.e., to what extent is the decay 
of a coherent state localized in the center of the regular 
island controlled by the RAT mechanism. Naturally, in 
cases wherein the regular islands are "resonance-free" the 
RAT mechanism is absent. A recent exa mple comes from 
the work of Feist et al. on their studies [187| on the con- 
ductance through nanowires with surface disorder which 
is controlled by dynamical tunneling in a mixed phase 
space. 

Thus based on the numerous systems studied there are 
clear indications that the nonlinear resonances are indeed 
important to understand the CAT mechanism. Neverthe- 
less certain observations reveal that one does not have a 
full understanding as of yet. For instance it is clear from 
Fig. [TJJ that the semiclassical theory underestimates the 
exact quantum splittings in the regimes of large Ti. This 
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could be due to the fact that the key nonlinear resonance, 
on which the estimate is based, is not clearly resolved. It 
could also be due to the fact that within the pendulum 
approximation the action dependence of the Fourier co- 
efficients are neglected. Alternatively it is possible that 
a different mechanism is operative at large fi and in par- 
ticular the Podolskiy-Narimanov theory might be rele- 
vant in this context. However, as ev ident from the re- 
sults reported in the recent work |l84l | by Mouchet et al., 
the Podolskiy-Narimanov estimate tends to overestimate 
in certain systems. In t his context note that Sheinman 
has recently derived 



|, using assumptions identical to 
those of Podolskiy and Narimanov, an expression for AE 
which differs from the result in Eq. (|32|) . Moreover the 
assumption of neglecting the action dependencies of the 
Fourier coefficients in Eq. (|55|) needs to be carefully stud- 
ied especially in cases wherein the nonlinear resonances 
are rather large. 

A more crucial issue pertains to the role played by 
the partial barriers like cantori at the regular-chaotic 
border. At present the RAT based mechanism still as- 
sumes the border to be devoid of such partial barriers 
and thus models the chaotic part by a random hermitian 
matrix. As briefly mentioned in the introduction the par- 
tial barriers do play an important role in the molecular 
systems as well. In this context the studies (l89l Il90j | 
by Brown and Wyatt on the driven Morse oscillator are 
highly relevant. The driven Morse system, a model for 
the dissociation dynamics of a diatomic molecule, clearly 
points to the cruc i al ro le of the cantori in the quan- 
tum dynamics 1 189L Il9dj . Interestingly the dissociation 
dynamics of the driven Morse system can be analysed 
from a dynamical tunneling perspective to understand 
the interplay of nonlinear resonances, chaos, and cantori 
in the underlying phase space. Investigations along these 
lines are in progress. Currently there is no theory which 
combines the RAT viewpoint with the effect of the hier- 
archical str uctur es located in the border regions of the 
phase space [191 1 . 



CONSEQUENCES FOR IVR 

Notwithstanding the issues raised in the previous sec- 
tions, the key point in the context of this review has 
to do with the relevance of the RAT and CAT phenom- 
ena to IVR. Heller in his stimulating review [l3( raised 
several questions regarding the manifestation of dynam- 
ical tunneling in molecular spectra. It is thus interesting 
to ask as to what extent has the recent progress con- 
tributed towards our understanding of quantum mech- 
anisms for IVR. This section therefore summarizes the 
recent progress from the molecular standpoint. 

The Hamiltonians used in this section are effective 
spectroscopic Hamiltonians of the form given in Eq. ([3]). 
The origin and utility of the effective Hamiltonians have 



already been discussed in the introduction. Three rea- 
sons are worth stating at this juncture. Firstly, the classi- 
cal limit of the spectroscopic Hamiltonians are very easily 
constructed. Secondly, the various nonlinear resonances 
are specified and hence it is possible to perform a careful 
analysis of the role of the resonances and chaos in dy- 
namical tunneling. A third reason has to with the fact 
that effective Hamiltonians of similar form arise in dif- 
ferent conte xts like electron transfer through molecular 
bridges [l92j |. coupled Bo se- Einstein con densates, and dis- 
crete quantum breathers [ 1 93l 1 1 94 1 1 951 ] . Interestingly the 
notion of local modes has close connections with the ex- 
istence of discrete breathers or intrinsic localized modes 
which have been, and conti nue to be, studied by a num- 
ber of researchers 1961 198l | . Wether the intrinsic local- 
ized modes can exist quantum mechanically is nothing 
but the issue of dynamical tunneling in such discrete 
systems. Therefore one anticipates that techniques of 
the previous sections should be useful in studying the 
breather-breather interactions as w ell. Indeed a recent 
work by Dorignac et al. invokes |l99l | the appropriate off- 
resonant states to couple the degenerate states. 



Local mode doublets 

To illustrate the RAT and CAT phenomena re- 
cent studies focus on the effective local mode 

Hamiltonian 12001 



H = H + g'vt? + + ^(^r + (49) 

with g' = g + X'(n\ + n 2 + 1) + A"n&, describing the 
highly excited vibrational states of a ABA triatomic. In 
the above the local mode stretches are denoted by (1,2) 
and the bending mode is denoted by b. As stated previ- 
ously such systems were one of the first ones to highlight 
the manifestation of dynamical tunneling in molecular 
spectra. Here, however, the bending mode is included 
and thus allows for several anharmonic resonances and a 
mixed phase space. The zeroth-order anharmonic Hamil- 
tonian is 

Ho = uj s {n\+ n 2 ) + uJbnb + x s {n\ + nf) + x\,n\ 

+ x sb n b (ni + n 2 ) + x ss n x n 2 (50) 

with rij = a/jdj being the occupation number of the j-th 

mode written in terms of the harmonic creation (aj) and 
the destruction (a,j) operators. The zeroth-order states 
\ni, n 2 , rib) are coupled by the anharmonic resonances of 
the form 



,(12) 







r(lb) 



r(26)^ 



v£p = (aln aj y + (atr( ai y 



(51) 



Thus VpxP couples states |n) and |n') with |nj — rii\ — q 



and In' 



p. Due to the structure of the Hamilto- 



nian the polyad number P = (m + n 2 ) + rib/ 2 is exactly 
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Table I: 


Local mode splitting 


i in cm 1 for H2O 


State" 




AE b ° 


(10 ± )0 


98.876 


100.26 


(20 ± )0 


48.278 


48.407 


(30 ± )0 


13.669 


12.919 


(40 ± )0 


2.661 


2.334 


(50 ± )0 


0.442 


0.514 


(60 ± )0 


0.105 


0.134 


(70 ± )0 


0.145 


0.013 



"Symmetrized lo cal m ode [nm labels 
'Taken from ref. l20ll 

c Calculated using the Baggott Hamiltonian. 



conserved and thus the system has effectively two degrees 
of freedom. 

The Hamiltonian has been proposed by Baggott [200j | 
to fit the vibrational states of H2O and D2O for specific 
choice of the parameters of H. A comparison of the local 
mode splittings predicted using t he Ba ggott Hamiltonian 
and the recent experimental data |20l| of Tennyson et al. 
is shown in the table HI The agreement is quite good 
considering the fact that the polyad P perhaps breaks 
down at higher levels of excitations. Another reason for 
the choice of the system had to do with the fact that the 
classical-quantum correspondence for the Hamiltonian in 
Eq. (|49l) was alr eady well established from several earlier 
works [202[ 12031 . l204l |. Additional motivation has to do 
with the possibility of studying the RAT and the CAT 
mechanisms in a controlled fashion by analysing the fol- 
lowing sub-Hamiltonians [206| : 



2V2 



(52) 
(53) 



The Hamiltonian Ha for smaller polyads exhibits near- 
integrable dynamics. The same is true for Hb but with 
an additional weak 2:2 resonance and hence it is pos- 
sible to study the influence of multiple nonlinear reso- 
nances. On the other hand the full Hamiltonian H has 
a mixed regular-chaotic phase space and thus on going 
from Ha — > Hb — > H a systematic study of the influ- 
ence of phase space structures on dynamical tunneling is 
possible. The classical limit can be obtained using the 
correspondence 



/ Ij exp(i0j 



Ij exp(- 



(54) 



and yields a nonlinear multiresonant Hamiltonian in 
terms of the action-angle variables of the zeroth-order 
Ho. For instance the resonant perturbation Eq. (f5l"j) has 
the classical limit 



V P -. q (I,0) 



(55) 
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Figure 12: (a) Local mode doublet splitting A,- in cm -1 with 
increasing excitation quantum r for polyad P = 6. The ex- 
act quantum (lines) are compared to the vibrational superex- 
change (symbols) at minimal order. Note that the agreement 
is excellent except for the r = 4 case for the full Hamiltonian 
(solid line). Inset shows the contribution from 130 differ- 
ent coupling chains towards obtaining As in case A. (b) A r 
estimated from the 1:1 nonlinear resonance alone (dashed) 
and the 1:1+2:2 resonances (dot-dashed) compared with the 
full case (solid line). Note that the Hamiltonian with 1:1 
and 1:1+2:2 are classically integrable. The three lower pan- 
els show the phase space (Ki/2,(j>i/4n) at increasing ener- 
gies corresponding to r — 2,3, and 4 for Eq. (|49|l . Signifi- 
cant amount of chaos and the 2:1 stretch-bend nonlinear res- 
onances set in for r = 4. 



with the obvious correspondence I — ► (n+ l/2)h. 

The zeroth-order states \n\ — r, ri2 = 0,n& = 2(P — 
r)) = |r, 0) and |0, r) are degenerate in the absence of 
the resonant couplings. However, in the presence of the 
couplings the degeneracy is broken resulting in the lo- 
cal mode doublets with a splitting A r . In other words 
if an initial state |r, 0) is prepared then on a time scale 
of T = 27r/A r the population is completely trasferred to 
the symmetric partner state |0, r). This corresponds to 
a dynamical tunneling process since th ere a re no ener- 
getic barriers to be blamed in this case |206| . Quantum 
mechanically for a given polyad P the splittings can be 
obtained easily by diagonalizing a N x N matrix with 
N = (P + 2)(P + 3)/2. The superexchange calculation 
in this multiresonant case is more involved than in the 
single resonance case. One important consequence of the 
multiresonant couplings is that the number of coupling 
chains connecting |r, 0) and |0,r) in t he state space for 
a given P is infinite. It turns out [205| that the minimal 
order estimate 
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provides a very good approximation to the exact quan- 
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turn splittings. The index /x labels all possible coupling 
chains T a i, c for a specific choice of a, b, and c satisfying 
the constraint a + 2b + 4c = 2r. The contribution to 
the overall splitting from a specific coupling chain, which 
has the form as in Eq. (23J), is denoted by A r (T^). The 
expression in Eq. (|56p is nothing but the lowest order 
of perturbation that would connect the two degenerate 
states involving a net change An = |n — n'| = 2r. The 
excellent agreement between the exact quantum and the 
superexchange calculation for P — 6 is shown in Fi g.H"2"R 
for the actual system and the various subsystems [206[. 
Note that in case of Ha a superexchange calculation of 
Ag involves twelfth-order perturbation with 130 coupling 
chains of varying signs, some as large as 10~ 4 , conspiring 
to yield a fairly accurate value. In Fig. [T2"b the splittings 
are calculated by retaining only the 1:1 resonance or the 
1:1+2:2 resonances. The A r monotonically decrease in 
both cases with increasing r i.e., OH-stretch excitations 
as observed earlier. However starting with r = 4 the inte- 
grable subsystem cannot account for the exact splittings. 
In particular the exact case exhibits nonmonotonic be- 
haviour which neither the integrable l imit case s nor the 
superexchange calculation can predict [205l l206j |. On the 
other hand the monotonic trend seen in Fig. [T2k for the 
near-integrable subsystem s is n ot true in general. For 
instance in an earlier work 205j strong fluctuations were 
seen in the near-integrable case for P = 8. It is well 
known that such fluctuations have to with the avoided 
crossings of the doublets with one or more states. The 
avoided crossings can occur in near-integrable (see Fig. [9] 
for example) as well as in the mixed phase space situa- 
tions. The phase space for the full Baggott Hamiltonian 
is shown in Fig. [12] with increasing energy. The su rface 
of sections are shown in the canonically coniugate [203j 
variables (ifi/2,0 x ) = ({h - 7 2 )/2, (61 - 2 )/2) at fixed 
sectioning angle <f> = 61+62 — 4:9b with <j> > 0. The normal 
mode regions appear around l\ ~ J 2 . The appearance of 
significant amount of chaos and the 2:1 stretch-bend res- 
onances in the phase space at higher energies suggests 
that the RAT and, to some extent CAT, mechanisms 
could explain the various results. 

As another example, in Fig. [13] the local mode split- 
tings computed using Eq. (f49|) for parameters appropri- 
ate to H2O and D2O and a higher polyad P = 8 are 
shown. It is clear that the splittings for a given local 
mode state are typically larger for D2O as compared to 
H2O. Note that this immediately establishes the dynam- 
ical nature of the tunneling since usual arguments of iso- 
topic substitution invoked for the potential barrier tun- 
neling cases would predict the opposite. This, of course, 
is not surprising since the dynamical tunneling involves 
transfer of vibrational quanta as opposed to the transfer 
of atoms. In addition the enhanced splittings for D2O 
is consistent with the fact that H2O is a better local 
mode system than D2O. The phase space for both sys- 
tems are mixed regular-chaotic and similar to that shown 
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Figure 13: Local mode splittings in cm -1 for the polyad P — 8 
for (A) H2O and (B) D2O using the appropriate parameters 
in the full Hamiltonian of Eq. (|49p . The open circles denote 
the exact quantum values whereas the filled circles represent 
a semiclassical calculation using an integrable approximation. 
Note that the integrable approximation decays monotonically 
in both cases and fails to account for the fluctuations. The 
phase space is mixed regular-chaotic for both the molecules 
but the fluctutations in A r are muted in D2O. 



in Fig. [12] Does the nonintegrability of the system dy- 
namics, and hence the chaos, have any impact on the 
local mode splittings? In order to address the question 
the splittings are computed for an integrable approxima- 
tion to the Baggott Hamiltonian. The splittings com- 
puted using the integrable approximation are compared 
to the exact splittings in Fig. [T3]for both the molecules. 
The integrable limit splittings, as expected, decay mono- 
tonically with increasing excitation. Although the com- 
parison looks favorable, except for the nonmonotonicity, 
for H2O it is clear from Fig. [13] that similar levels of 
agreement are not seen for D 2 0. More importantly, the 
fluctuations in A r are completely missed and thus sug- 
gest that the exact splittings are determined by a subtle 
interplay of the nonlinear resonances and chaos. 

To illustrate RAT consider the mixing between the 
zeroth-order states |2,0) and |0, 2) due to the stretch- 
bend resonance in Ha- The states are not coupled di- 
rectly by the resonances as can be inferred by the state 
space shown in Fig.H"4"a. The reso nanc e zones, computed 
using the Chirikov approximation [207J, are quite far from 
the states. Indeed Fig. [Mb shows that the classical dy- 
namics exhibits trapping, however the quantum survival 
probability of |2, 0) shows complete mixing with |0, 2) on 
a time scale of about 0.1 nanoseconds. As before off- 
resonant states participate in vibrational superexchange, 
clear from Fig. 1141 1. resulting in the observed mixing. 
There are now several coupling chains at minimal order 
and the perturbative calculation yields A 2 ~ 0.18 cm -1 
in excellent agreement with the exact value (cf. Fig. [T2"]) . 
What is the role of the stretch-bend resonances? This 



25 




as 1 20 



2 4 6 0.1 0.2 



CM 




.02 



Figure 14: (a) State space (111,11,2) for P = 6 with the loca- 
tion of the 2:1 stretch-bend nonlinear resonance zones (dashed 
regions) shown for reference, (b) Quantum survival proba- 
bility for the state |2, 0) (circles) computed using Ha shows 
two-level like mixing with the symmetric partner |0, 2) de- 
spite trapping of classical trajectory on the same timescale 
(inset). Thus the quantum energy transfer is due to dynam- 
ical tunneling, (c) RAT mechanism due to the 1:1 nonlinear 
resonance (shaded) in the phase space induced by the remote 
2:1s explains the energy flow, (d) Overlap intensities versus 
6E^ = (E a — E^) (in units of mean level spacing) shows in- 
volvement of off-resonant states. 



becomes clear on looking at the phase space for Ha at 
E = E 2 8 shown in Fig. WSz . One immediately notices 
a nonlinear resonance located between the phase space 
regions corresponding to the two states. The resonance 
however is not the 2:1 stretch-bend resonance but a weak 
1:1 resonance between the local stretch es in duced by the 
stretch-bend couplings. It ca n be shown |206l ] using classi- 
cal perturbation theorv |l80l ] that the induced 1:1 is gov- 
erned by a pendulum Hamiltonian of the form 



H(Kt,<f>i;K 2 ,N) 
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where the parameters 57 s , a ss , and a\ 2 are related to the 
original para meters of Eq. (|49|) . The induced coupling 
strength[208( is related to the width of the resonance is- 
land seen in Fig. [T4h . Thus having identified the key 
resonance the splitting between a general local mode pair 
can be estimated as 
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(59) 

= j 2 /2Afn being the unperturbed energy from 
Using parameters appropriate for H2O in 
and the expressions in Eq. (|58c() one obtains 
;' 1.32 x 1(T 2 and g md « 3.93 cm" 1 . Thus 
0.20 cm -1 which is in agreement with the ex- 



act and superexchange values. This clearly supports the 
RAT mechanism for the quantum mixing seen in Fig. 114b . 
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Figure 15: (a) Splittings in cm -1 for P = 6 overtone states for 
Ha- The exact result (solid line), superexchange calculation 
(open triangles), and predictions based on the induced 1:1 
stretch-stretch resonance alone (filed triangles, dashed line) 
are compared, (b) Phase space for Ha ~ ^6 emphasizes 
the need for invoking multiple resonances (arrows). RAT 
mechanism using the multiple resonances leads to very good 
agreement shown in (a) as shaded cricles. 



with K 2 = h+h (1:1 subpolyad) and N = (I 1 +I 2 )+I b /2 
(classical analog of the polyad P) being the approxi- 
mately and exactly conserved quantities. The parameters 
appearing in the pendulum Hamiltonian can be expressed 
in terms of the parameters of the original Hamiltonian 



Although the classical dynamics of Ha is near- 
integrable the phase space can be very rich in terms 
of several resonances. The RAT mechanism depends on 
identifying the crucial resonances. In Fig. [T^the mono- 
tonic decrease of A r with r for the system described by 
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Ha and Hb might suggest that the induced resonance 
alone controls the splittings. However Fig. [T5b shows 
that calculating the splittings using Eq. (|5T|) yields very 
poor results for r — 5,6. The reason becomes clear on 
inspecting the phase space and Fig. [T5b shows the case of 
r = 6 as an example. The Husimi distribution of the state 
|6, 0, 0) would be localized about K\ ss 6 and thus the 2:1 
stretch-bend and the 1:1 induced resonance islands seper- 
ate the state from its symmetric partner |0, 6, 0) localized 
around K\ ps — 6. Several other resonances are also visi- 
ble at this energy in Fig. 115b but they are of higher order 
and perhaps not relevant for large h va lues. Moreover 
the Husimi representations indicate [206j that the state 
|5,0, 2) is localized inside the 2:1 island (ground state) 
and the state |4, 0, 4) corresponds to the separatrix state. 
Thus the 2:1 resonance must be coupling the |6, 0, 0) and 
|4, 0, 4) states. As a consequence the RAT mechanism 
can be schematically written as 



|6,0) ^ |4,0) m |0,4) 2$ | ,6) 
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with f3 e ff being the effective coupling across the 2:1 is- 
lands and <7i n d being the effective 1:1 coupling in Eq. (|58c|) 
for r = 6. This is illustrated in Fig. [T5b on the phase 
space as well. Again approximating the 2:1 islands by 
pendula it is possible [206J to determine the /3 e ff and the 
computed splittings agree rather well with the exact re- 
sults as shown in Fig. [T5a . Thus Fig. [15] demonstrates 
the RAT mechanism in that the mixing between two re- 
mote (both in state space and the phase space) states 
|6, 0, 0) and |0, 6, 0) is characterized by the nonlinear reso- 
nances. Further confirmation of the multiresonant mech- 
anism comes from the fact that values of Aq for Hb and 
the full Hamiltonian also compare well with the exact re- 
sults. The subtelty of the mechanism lies in the fact that 
all the nonilinear resonances in Fig. 115b arise from the 
two 2:1 resonances and at smaller values of h one expects 
the higher order resonances to come into play. 

The role of RAT in explaining the local mode doublet 
splittings is obvious from the preceeding discussion and 
examples. Model studies reveal that even very small, 
on the scale of h, nonlinear resonan ces m ediate dynam- 
ical tunneling in molecular systems (206] . However the 
role of CAT is far from being established and some of 
the issues involved were highlighted before in this re- 
view. In essence, to implic ate CAT a tleast a third chaotic 
state has to be identified 147 , Il48| . Moreover the sys- 



tem should be in the deep semiclassical limit i.e., small 
enough H so that the quantum dynamics can sense the 
presence of chaos in phase space. These are fairly severe 
constraints and experimentally there are not many sys- 
tems which respect the constraints. Intutively speaking, 
highly excited energy regions of moderate sized molecules 
might be ideal systems. One possibility is that the sig- 
natures of CAT could be present in the eigenstates of 
highly excited molecules and could potentially compli- 
cate the state assignments. One such example can be 
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Figure 16: Variations in the splittings for the ground state of 
the 2:1 island with coupling strength f3. (a) Effect of breaking 
the polyad by adding weak 3:1 resonance to H with strength 
5 — 0.1 cm -1 (circles) contrasted with the polyad conserved 
(P — 16) case (thick line). The phase space (with axes as 
in Fig. I12p is shown in the latter case for select values of /?. 
(Inset) Fluctuations in the splitting for P — 8 case (thick 
line) compared to that with Ti reduced by half (circles). All 
quantities are in units of cm -1 . 



found in the work [202} by Keshavamurthy and Ezra on 
the dynamical assignments of the highly excited states 
of H2O for P — 8. Another example, perhaps, can be 
found in the work |209l ] by Jacobson et al. wherein the 
highly excited (~ 15000 cm -1 above the ground state) 
vibrational states of acetylene are dynamically assigned. 
In their study Jacobson et al. found that approximate 
assignments could be done in an energy range which had 
highly chaotic phase space whereas many of the eigen- 
states could not be assigned in another energy range de- 
spite the existence of large regular regions in the phase 
space. To some extent an indirect role of CAT is evi- 
dent from Fig. [T2"b wherein the dominant integrable res- 
onances cannot account for the exact splittings. Some 
more hints originate from Fig. 1161 where the fluctuations 
in the splittings associated with the ground state of the 
2:1 island are shown in two different situations. The first 
of these investigates the effect of breaking the polyad con- 
stant P by adding weak 3:1 stretch-bend resonances to 
the Hamiltonian in Eq. (|4"5|) . Thus the three degrees of 
freedom Hamiltonian is 



H 3d = H + S(v£ b) 



(61) 



with 5 — 0.1 cm 1 which is three orders of magnitude 
smaller than the mean level spacing of H. A consequence 
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of 6 7^ is that the effective density of states is increased 
providing many more states with different P to interact 
with the state of interest. In Fig. [THlthe effect of such a 
small perturbation is shown for P = 16 and clearly there 
are far more fluctuations in the mixed regime. However 
away from the fluctuations the inverse participation ratio 
of the state is nearly one in the basis of the two degrees of 
freedom Hamiltonian H which suggests minor influence 
of states with P ^ 16. Despite this and the weak cou- 
pling strength S there are significant differences between 
the computed splittings for H^d and H. Presumably this 
has to do with the subtle changes in the nature of the 
chaotic states upon adding the 3:1 resonances but further 
studies are required in order to gain a better understand- 
ing. The results of a second system shown in the inset to 
Fig. [TH] reveal the increased flutuations in the splittings 
for the state with P = 8 on decreasing the value of % 
by half. Interestingly in this case the average splitting 
seems to decrease with increasing (3 and thus increasing 
stochasticity in the phase space. 

Most of the detailed phase space studies have focused 
on the role played by classical phase space structures for 
vibrational modes alone. The involvement of rotations, 
torsions, and other large amplitude modes in m olec- 
ular systems and their implications for IVR [2ld | and 
dynamical tunneling have received relatively less atten- 
tion. Early pioneering studies 0] by Harter and Patterson 
showed that the splitting of the doubly degenerate K- 
levels of a symmetric top on the introduction of asymme- 
try could be associated with tunneling across a separatrix 
on the so called rotational energy surface. Would the cou- 
pling of vibrations to these other modes enhance or sup- 
press the splittings? Harter has extended the rotational 
energy surfaces idea to describe and interpret the dynam- 
ics and spectra for a class of anharmonic coriolis cou- 
pled vibrational modes 118 ] . Lehmann studied [21 ll |212| | 
the coupling of rotations to the local mode vibrations 
in XH„ systems. In p arti cular Lehmann, and indepen- 
dently Child and Zhu [213| . considered the modification 
of the local mode dynamical tunneling due to the in- 
teractions with rotations. Depending on the local mode 
tunneling time the rotation motion would reflect the ex- 
tent of dynamical symmetry of the local mode state. In- 
terestingly Lehmann's analysis reveals that it is possible 
for the rotation of a molecule to significantly redu ce th e 
tunneling rate i.e., decrease the doublet splitting 211 1 . 
Such rotational quenching of the local mode tunneling is 
due to the fact that in addition to the transfer of vibra- 
tional quanta it is also necessary to reorient the angular 
momentum in the body fixed frame. Semiclassical in- 
sights arise by looking at the motion on the rotational 
energy surfaces. It is not clear as to how much of a role 
does the classical chaos plays in the studies of Harter 
and Lehmann. However, the local mode splitting is de- 
termined by RAT and CAT mechanisms and hence they 
also determine the timescales relative to rotations. It 



remains to be seen if the different mechanisms are re- 
flected in the nature of the rovibration al eig enstates. A 
possible signature of CAT was suggested [2 161 ] by Ortigoso 
nearly a decade ago. Ortigoso found anomalously large 
iT-splittings for some of the asymmetry doublets occur- 
ing in the rotation-torsion energy levels of acetaldehyde. 
Further studies involving coupling with vibrations have 
not been undertaken although, based on time scale sep- 
arations one can argue that the observed enhancements 
might lead to increased vibration-rotation couplings with 
important consequences for IVR. 



State mixing in multidimensions 

Every single example in this section up until now has 
involved systems with discrete symmetry and two degrees 
of freedom. The requirement of discrete symmetries for 
dynamical tunneling is not very restrictive. As discussed 
before in this review this aspect has b een emphasized in 
the early work as well. Tomsovic[218j has also discussed 
the possibility of observing CAT in the absence of reflec- 
tion symmetries. On the other hand studies on dynami- 
cal tunneling in systems with large degrees of freedo m ar e 
crucial. Given that there are important differences [H^], 
briefly mentioned in the introduction, between the clas- 
sical dynamics in two and more than two degrees of free- 
dom it is natural to ask if the RAT and CAT mechanisms 
hold in general. Especially from the IVR perspective the 
crucial question is wether dynamical tunneling can pro- 
vide a route for mixing between near-degenerate states 
and hence energy flow between phase space regions sup- 
porting qualitatively different types of motion. Thus is it 
possible that an excited CH-stretch can, over long times, 
evolve into a CC-stretch despite the lack of an y obvi- 
ous coupling? Heller and Davis coni ect ured 1 1 1 4l ] that it 
should be possible and the agent would be dynamical 
tunneling. However support for the conjecture from a 
phase space viewpoint have not been forthcoming. 

A first step towards such goals was recently taken [22C| 
by studying the model Hamiltonian 



H = Hn 



(a{)^(a 2 )^(a 3 )^(a^ + h. 



describing four coupled modes j 
zeroth-order part 



(62) 

1,2,3,4 with the 



H = ^(LUjUj 



) + x ^njn k (63) 
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The zeroth-order states |n), eigenstates of Ho with 
eigenvalues E^, get coupled and hence mixed due 
to the anharmonic resonances chracterized by m r = 
(av, — j3 r , — Y r , — S r ) with strengths K mr . In order to 
model a system without any symmetr ies t he parame- 
ters of Hq were taken from the work 22 1| of Beil et 



28 



0.75 




0.25 



3400 



3600 



3800 



E°/AE 
n 



Figure 17: (A) Dilution factors a of the zeroth-order states 
jn) with H \n) = E°\n) versus energy for the Hamiltonian in 
Eq. (|62p . The polyad value is P = 8 and couplings 0.5-ftTmi = 
Km 2 = Km 3 ~ 0.2AE with mean level spacing AE « 4.4 
cm -1 . (B) Histogram showing the distribution of a in (A). 
Note that a large number of states have unit dilutions but 
there are also a number of states with a < 1 indicating mixing 
of the zeroth-order states. Is the mixing due to dynamical 
tunneling? In (C) an expanded view of the states around 
E w 3542. 5AE in (A) is shown. Four states of interest are 
highlighted by squares. 
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Figure 18: (Colour online) Quantum (top panel) and classi- 
cal (bottom panel) survival probabilities of the four zeroth- 
order states shown in Fig. I11C . Parameters chosen are as in 
Fig. 1111 Time r is measured in units of the Heisenberg time 
th = (2-kcAE)^ 1 . State \a) (black) shows long time trap- 
ping in both cases. States |6) (red), |c) (green), and \d) (blue) 
mix amongst each other quantum meachanically in contrast 
to their classical behaviour. The quantum cross probabili- 
ties |(n'|n(t))| 2 ar e als o shown with consistent color coding. 
Adapted from ref. |220| . 



al. corresponding to the four high frequency modes 
of the molecule CDBrClF. Three perturbations mi = 
(1,-2, 0,0), m 2 = (1,-1,-1,0), and m 3 = (1,-1,0,-1) 
with very weak strengths K mr /AE = k mr < 1, rela- 
tive to the mean level spacing AE, are selected. Note 
that the coupling structure of the Hamiltonian implies 
the existence of a good polyad quantum number P = 
rii + (ri2 + na + n^)j1. Thus the model Hamiltonian has 
effectively three degrees of freedom. 

A few words about the choice of the Hamiltonian is 
perhaps appropriate at this stage. Firstly, the form of 
Hq is quite generic and t he parameters were chosen from 
an experimental fit 22l\ in order to test the ideas on 
a real molecular system. Secondly the Hamiltonian in 
Eq. (|rj2")) can be cons idere d to be in the intrinsic res- 
onance representation 222j . In other words the Hamil- 
tonian models the weak residual couplings between the 
KAM tori corresponding to Hq. Finally the choice of 
the resonant couplings is also arbitrary; a different set of 
couplings would generate mixings among different zeroth- 
order states. 

The relevant issue here has to do with the fate 
of zeroth-order states |n), |n'), . . . , which are near- 
degenerate E® ss in the presence of such weak 
perturbations. In Fig. [TTR the dilution factors for the 
various zeroth-order states are shown and it is clear that 
many zeroth-order states are mixed. In particular within 



one mean level spacing several states have varying levels 
of mixing and a typical example is shown in Fig. I17C . 
Is the mixing observed in Fig. I17C due to dynamical 
tunneling? In order to check four states are picked and 
correspond to \a) = |0, 11, 1, 4), \b) = |0, 11, 2, 3), |c) = 
|0, 12, 2, 2), and \d) = |0, 13, 1, 2). The choice of the en- 
ergy range and the states is fairly arbitrary as evident 
from Fig. [TTJ First it can be confirmed by evaluating the 
quantum survival probability P„q and its classical ana- 
log P n c, shown in Fig. [TH1 that the mixing of the states 
of interest is indeed a quantum process i.e., classically 
forbidden. It is also apparent from Fig. [TS] that states 
| b), |c), and \d) mix amongst each other over long time 
scales. State \a) on the other hand shows trapping both 
classically and quantum mechanically. Explanation for 
the observed mixings in Fig. [T7] can be given in terms 
of the vibrational superexchange mechanism. Indeed the 
IVR 'fractionation' pattern, shown in Fig. Q33 indicates 
a a clump of virtual or off-resonant states about 60Ai5 
away from the corresponding 'bright' states. Such frac- 
tionation patterns are similar to those observed in the 
experiments [9^] of Callegari et al. for example. The off- 
resonant states have a ~ 1 and hence do not mix signifi- 
cantly. However it is important to note that the fraction- 
ation pattern looks the same in case of every state and 
hence there must be subtle phase cancellation effect to 
explain the trapping exhibited by \a). At the same time 
Fig. [TH] also shows that the vastly different IVR from the 
four states cannot be obviously related to any avoided 
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crossing or multistate interactions. 
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Figure 19: (Colour online) The top panel shows variation 
of the eigenvalues E = E a /AE with the coupling parame- 
ter fe mi with the other couplings fixed as in Fig. 1111 Eigen- 
states having the largest contribution from the specific zeroth- 
order states are indicated. Bottom panel shows the 'frac- 
tionation' pattern i.e., overlap intensities p na versus Ae n = 
(-En — E a )/AE on log scale. The cluster of states around 
Ae n « 60 are the off-resonant states that mediate the mixings 
seen in Fig. [17] and Fig. [18] via a vibrational superexchange 
mechanism. Adapted from ref. |220| . 



In light of the conjecture of Davis and Heller [114j and 
the recent progress in our understanding of the mech- 
anistic aspects of dynamical tunneling it is but natu- 
ral to associate the near-degenerate state mixings with 
RAT. However how does one identify the specific reso- 
nances at work in this case? Given that the system has 
three degrees of freedom it is not possible to visualize the 
Poincare surface of section at H E. In the pr esent case 
the cruicial object to analyse is the Arnol'd web |l8f| i.e., 
the network of nonlinear resonances and the location of 
the zeroth-order states therein. Fortunately the classical 
limit Hamiltonian 



H{1,6) = H (I) + 2eJ2 K ™ r \/ I ? rI 2 rI 3 r 4 r cos(m r 6) 

r 

= H (I)+eJ2 K mJAl)cos(m r -9) (64) 

r 

from Eq. (|62p is easily obtained, as explained earlier, in 
terms of the action-angle variables (1,9) of Hq. In the 
above equation the variable e is introduced for conve- 
nience during a perturbation analysis of the Hamiltonian. 
The 'static' Arnol'd web at E w E and fixed polyad 
P c = I\ + (1% + h + h)/2, classical analog of the quan- 
tum P , can then be constructed via the intersection of 



the various resonance planes m r ■ dHoiVf/dl = with 
the energy shell Tto(T) *=* E. The static web involves all 
the nonlinear resonances restricted to, say, some maxi- 
mum order \a r \ + \/3 r \ + \j r \ + \5 r \ < M. The reason for 
calling such a construction as static has to do with the 
fact that being based on Hq it is possible that many of 
the resonances do not have any dynamical consequence. 
Thus, although the static web provides useful informa- 
tion on the location of the various nonlinear resonances 
on the energy shell it is nevertheless critical to determine 
the 'dynamical' Arnol'd web since one needs to know as 
to what part of the static web is actually relevant to 
the dynamics. Further disc ussio ns on this point can be 
found in the paper by Laskar[223] wherein time-frequency 
analysis is used to study transport in the four dimen- 
sional standard map. In order to construct the dynam- 
ical Arnol'd web it is necessary to be able to determine 
the system fr eque ncies as a function of time. Several 
techniques [53, 223 1 have been suggested in the literature 
for performing time-frequency analysis of dynamical sys- 
tems. An early example comes from the work |224j of Mil- 
czewski, Diercksen, and Uzer wherein the Arnol'd web for 
the Hydrogen atom in crossed electric and magnetic fields 
has been computed. A critical review of the various tech- 
niques is clearly outside the scope of this work and hence, 
without going into detailed discussion of the a dvantages 
and disadvantages, the wavelet based approach 225l 226| 
developed by Wiggins and coworkers is utilized for con- 
structing the web. There is, however, ample evidence 
that the wavelet based local frequency analysis is ideally 
suited for this purpose and therefore a brief description 
of the method follows. Classical trajectories with initial 
conditions satisfying H(I, 9) w E are generated and the 
nonlinear frequencies Q k (t) are computed by perform- 
ing the continuous wavelet transform of the time series 
z k (t) = y/2I k (t) exp(t0 fc (t)): 



L g Zk(a,b) = a 1/2 



Zk(t)g* 



dt 



(65) 



with a > and real b. Various choices can be made for 
the mother wavelet g(t) and in this work it is chosen to 
be the Mor let- Grossman function 



l(t) 



V2 



1 f n x t 2 

exp I 2m Xt — — j 



(66) 



and the parameter values A = 1, a = 2. Note that 
Eq. (|6"5|) yields the frequency content of z k within a time 
window around t = b. In many instances one is in- 
terested in the dominant frequency and hence the re- 
quired local frequency is extracted by determining the 
scale (a, inversely proportional to frequency) which max- 
imizes the modulus of the wavelet transform [225j i.e., 
n k (t = b) = max a \L g z k (a,b)\. The trajectories are fol- 
lowed in the frequency ratio space (Oi / O3, Ox/04) since 
the energy shell, resonance zones, and the location of 
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the zeroth-order states can be projected onto the space 
of two independent frequency ratios. Such 'tune' spaces 
have been constructed and analysed before, for instance, 
by Martens, Davis, and Ezra in their seminal work[53j on 
IVR in the OCS molecule. A density plot is then created 
by recording the total number of visits by the trajectories 
to a given region of the ratio space. Quite naturally the 
density plot highlights the dynamically significant por- 
tions of the Arnol'd web at a given energy and polyad. 
The computation of such a dynamical web is shown in 
Fig. [201 at an energy E corresponding to the zeroth-order 
states of interest. One immediately observes that apart 
from highlighting the primary resonances Fig. [20] shows 
the existence of higher order induced resonances. More 
importantly the zeroth-order states are located far away 
from the primary resonances and thus the state mixings 
cannot be ascribed to the direct couplings in the Hamil- 
tonian. However, the states are located in proximity to 
the junction formed by the weak induced resonances de- 
noted as mil, m^, and m^. The nature of these induced 
resonances make it very clear that the state mixings are 
due to RAT. 




2.1 2.2 2.3 2.4 



"1 /U 3 

Figure 20: (Colour online) The dynamical Arnol'd web at 
E = E generated by propagating 25000 classical trajecto- 
ries for r « 40. The primary resonances mi (red), m2 
(black), and 1113 (magenta), as predicted by Ho are super- 
imposed for comparison. The zeroth-order states of interest 
however are located far away from the primary resonances 
and close to the junction of three weak induced resonances 
mil = (0, 0, 1, -1), m l2 = (0, 1, -1, 0), and m« = (0, 1, 0,-1) 
indicated by arrows. The induced resonances lead to the ob- 
served state mixings in Fig. 1 171 and Fig. [19] via the mechanism 
of RAT. The data for this fig ure w ere generated by A. Sem- 
parithi. Reproduced from ref. |220| . 



How can one be sure that it is indeed the RAT mecha- 
nism that is at work here? The simplest answer is to try 
and int erfer e with the RAT mechanism by removing or 
altering [220| the relevant resonances as seen in Fig. [501 



This was precisely what was done in the recent work [220] 
and classical perturbation theory is again used to remove 
the primary resonances in Eq. ([64]) to O(e). Since the 
methodology is well known only a brief description of 
the perturbative analysis follows. The first step involves 
reduction of the Hamiltonian in Eq. (|64[) . noting the ex- 
actly conserved polyad P c , to an effective three degrees of 
freedom system. This is achieved by performing a canon- 
ical transformation (I, 9) — > (J, if); P c ) via the generating 
function 

F(J,^:P c ) = Y,(^r-e)J r + e 1 P c (67) 

r 

The new variables can be obtained in terms of the old 
variables by using the generating function properties ip = 
djF and I = dgF. The reduced Hamiltonian takes the 
form 

H(3, 4>; P c ) = H {3; P c ) + e ]T K mr g r (3; P c ) cob W 

r 

(68) 

with g r (3; P c ) being determined from the functions f r (T). 
The resonant angles are related to the original angle vari- 
ables by ip r = m r • 9. The dynamical web in Fig. [20] 
clearly shows that the zeroth-order states of interest are 
far away from the primary resonances. Thus the three 
primary resonances are removed to 0(e) by making use 
of the generating function 

G = ^ tp r J r + e d r (3) sini/v (69) 

r r 

where the unknown functions d r {3) are to be chosen such 
that the reduced Hamiltonian does not contain the pri- 
mary resonances to O(e). Using the standard generating 
function relations one obtains 

J r = J r + ed r (3) cosi/v (70) 
ij> r — tp r + e d rr r(3) sin-i/v (71) 

r' 

with r = 1,2,3 and d rr >(3) = dd r (3)/dJ r i . Writing the 
reduced Hamiltonian in terms of the variables (J, xp) and 
using the identities involving the Bessel functions J n 

cos(eet sin b) = Mea) + 2^J 2l {ea)cos{2lb) (72) 

i>i 

sin(easin6) = 2^J 2 i+i sin((2Z + 1)6) (73) 
z>o 

one determines the choice for the unknown functions to 
be 

dr{ j ) = - Kmr 9r&3± (74 ) 

Consequently, at 0(e 2 ) the effective Hamiltonian 

H(3, 4>; P c ) w H (3; P c ) + e 2 H 2 (3, -0; P e ) (75) 
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containing the induced resonances is obtained. 

Following the procedure for RAT the induced reso- 
nances are approximated by pendulums. For example 
in the case of one obtains 



n 



(24) _ (if 24 - K^f 



eff 



2M- 



21 



+ 2V24Cos(202 4 ) (76) 



with K24 ~ 2/4 and 204 ~ (#2 — 6* 4 ). The resonance 
center is and the effective coupling can now be ex- 
pressed in terms of the conserved quantities l3,P c , and 
P24 = h + Ii- Note that P24 is the appropriate polyad for 
the induced resonance 172:^4=1:1. Thus it is possible to 
explicitly identify (cf. Eq. (|76|)) the barrier for dynamical 
tunneling in this three degrees of freedom case. Observe 
that the barrier is parametrized by an exact constant of 
the motion P c and two approximate constants of the mo- 
tion I 3 , and P 2 4- The effective couplings can be trans- 
lated back to effective quantum strengths A m V m /2 
quite easily. The perturbative analysis, for parameters 
relevant to Fig. [TSl yields 



Amu < l A 



m i3 I 



0.07K n 



(77) 



and thus the induced resonances are more than an or- 
der of magnitude smaller in strength as compared to the 
primary resonances. Based on the RAT theory it is pos- 
sible to provide an explanation for the drastically differ- 
ent behaviour of state \a) when compared to the other 
three states as seen in Fig. Q1JJ The state \a) is not sym- 
metrically located, refering to the arguments following 
Eq. I|38p. with respect to |6) and thus, combined with the 
smallness of A mil , does not show significant mixing. 

Given that the key resonances and their strengths are 
known (cf. Eq. ([77|) ) is it possible to interfere with the 
mixing between the states? In other words one is ex- 
ploring the possibility of controlling the dynamical tun- 
neling, a quantum phenomenon, by modifying the local 
phase space structures. Consider modifying the quantum 
Hamiltonian (based on classical information!) ofEq. (|62|) 
as follows: 

H' = H+\X m J(ala 3 + H.c.) + \\ m J(ala 4 + H.c.) (78) 

In the above terms have been added to counter the in- 
duced resonances. Due to the weakness of the induced 
couplings quantities like mean level spacings, eigenvalue 
variations shown in Fig.ll9( and spectral intensities show 
very little change as compared to the original system. 
However Fig. I2T1 shows that the survival probabilities of 
all the states exhibit long time trapping. The dilution 
factors shown in Fig. [21] also indicate almost no mix- 
ing between the near-degenerate states. Thus dynam- 
ical tunneling has been essentially shut down for these 
states. The fractionation pattern for the states with the 
modified Hamiltonian is quite similar to the ones seen 
in Fig. []j|J One again observes a clump of off-resonant 



states around the same region. Surely the vibrational su- 
perexchange calculation would now predict the absence 
of mixing due to subtle cancellations but it is clear that a 
more transparent explanation comes from the RAT mech- 
anism. More significantly Fig.l2"T1shows that other nearby 
zcorth-order states are unaffected by the counter reso- 
nant terms. Thus this is an example of local control of 
dynamical tunneling or cquivalcntly IVR. 




3542 3542.5 3543 
E>E 

Figure 21: (Colour online) The top panel shows the sur- 
vival probabilities for the states of interest using the modi- 
fied Hamiltonian in Eq. (|78p . Note that the counter resonant 
terms have essentially shut off the dynamical tunneling prov- 
ing the RAT mechanism. The bottom left figure shows the 
change in the dilution factors. Solid symbols refer to the 
a calculated from the original Hamiltonian in Eq. (|62[) (cf. 
Fig. \17C ) and the open squares are the a calculated using 
Eq. (|78p . Only the states influenced by the induced reso- 
nances in Fig. [20] have their a — *• 1 whereas nearby states 
show very little change. The four panels on the right show 
the overlap intensity for the modified Hamiltonian. Notice 
the clump of states near m 60A_E as in Fig. 1191 



The model Hamiltonian results in this section corre- 
spond to the near-integrable limit and thus RAT ac- 
counts for the state mixings. It would be interesting to 
study the mixed phase space limit for three degrees of 
freedom systems where classical transport mechanisms 
can compete as well. This however requires a careful 
study by varying the effective h of the system in order 
to distinguish between the classical and quantum mech- 
anisms. Note that such zeroth-order state mixing due to 
weak couplings occurs in other model systems as well. 
In particular the present analysis and insight suggests 
that the mixing o f edg e and interior zeroth-order states 
in cvanoacetvlene [106j must be due to the RAT mecha- 
nism. Two more examples in Fig.[22]highlight the mixing 
between near-degenerate states due to very weak cou- 
plings. The case in Fig. I2"2"a corresponds to the mixings 
induced by the weak 3:1 resonance in the case of the 
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Figure 22: Additional examples for state mixings due to weak 
couplings, (a) Interpolyad mixing between the eigenstates of 
the Baggott Hamiltonian due to the polyad breaking Hamilto- 
nian in Eq. (|61|l with 8 — 0.1 cm - . The inverse participation 
ratios of the various eigenstates of H^d in the Baggott eigenba- 
sis are shown versus the eigenvalues E in cm -1 , (b) Dilution 
factors of the zeroth-order states for a Hamiltonian of the 
form as in Eq. (|62|l but with a different Ho- The polyad cho- 
sen is P = 16 and the three resonant couplings have strengths 
0.5K mi = -Km 2 = Km 3 ~ 0.05 in units of the mean level spac- 
ing, AE w 1.8 cm -1 . The zeroth-order energies E are also in 
units of AE. 



model Hamiltonian of Eq. (|6ip with 5 = 0.1 cm The 
inverse participation ratios of the eigenstates |cv3d) of H^d 



L (3d) 



(79) 



are calculated in the Baggott eigenbasis {|ap)} spanning 
several polyads centered about P = 16. By construc- 
tion the participation ratio of every eigenstate would be 
close to unity if the small 8 value did not induce any 
interpolyad mixing. However it is clear from Fig. I2"2"k 
that several eigenstates do get mixed. The other ex- 
ample shown in Fig. 122b pertains to a case where the 
Hamiltonian is of the same form as in Eq. (|62[) but with 
the pararnetrs of Hq that corresponds to the CF3CHFI 
molecule [227| . The dilution factors of the various zeroth- 
order states in the polyad P = 16 with extremely small 
couplings exhibit mixing. In both the examples it would 
be interesting to ascertain the percentage of states that 
are mixed due to CAT/RAT. 



SUMMARY 

The central message of this review is that both clas- 
sical and quantum routes to IVR are controlled by the 
network of nonlinear resonances. A simple sketch of 
IVR in state space via classical diffusion and dynamical 




Figure 23: (Colour online) A sketch of the resonance web 
in three degrees of freedom in the (Ii, 12,13) action (state) 
space. The energy surface Ho = E is intersected by several 
resonances which form the web and control the IVR out of 
a state located on the energy surface. The strength of the 
resonances is indicated by the thickness and shaded regions 
highlight some of the resonance junctions. Mixing i.e., en- 
ergy flow between the state on the top (red circle) and a far 
away state (blue square) is due to RAT, CAT (a possible tun- 
neling sequence is shown by red arrows), classical across and 
along (thick arrows) transport. The competition between the 
tunneling and classical transport mechanism depends on the 
effective Ti of the system. Barriers to classical transport do 
exist and are expected to play an important role[53l] but they 
have not been indicated in this simple sketch. 



tunneling is shown in Fig. [23] and should be contrasted 
with the sketch shown in Fig. [3] The vibrational su- 
perexchange picture of Fig. [3] explaining IVR through 
the participation of virtual off-resonant states is equiva- 
lent to the picture shown in Fig. [23] which explains IVR 
in terms of the Arnol'd web and states on the energy 
shell. Furthermore, as sketched in Fig. [23] a given zeroth- 
order state is always close to some nonlinear resonance 
or the other. Thus it becomes important to know the 
dynamically relevant regions of the web. Classically it 
is possible to drift along the resonance lines and come 
close to the junction where several resonances meet. At 
a junction the actions can move along an entirely dif- 
ferent resonance line. In this f a shio n, loosely referred 
to as Arnol'd diffusion [2281. [229l [230l |. the classical dy- 



namics can lead to energy flow between a given point 
in Fig. [23] and another distant point. Arnol'd diffusion, 
however, is exponentially slow and its existenc e require s 
several constraints to be satisfied by the system 23lll232 |. 
Therefore dynamical tunneling between two far removed 
points shown in Fig. [23] is expected to be more proba- 
ble and faster than the classical Arnol'd diffusion. On 
the other hand, for generic Hamiltonian systems, a spe- 
cific primary resonance will be intersected by infinitely 
many weaker (higher order than the primary) resonances 
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but finitely many stronger resonances. The stronger res- 
onances can lead to a characteristic cross-resonance dif- 
fusion near the resonance junction. The cross-resonance 
diffusion, in contrast to the Arnol'd diffusion, t akes p lace 
on time scales much shorter than exponential [2331 ] . In 
Fig. [221 one can imagine a scenario wherein a sequence 
of dynamical tunneling events finally lands up close to a 
resonance junction. Would the classical erratic motion 
near the junction interfere with the dynamical tunnel- 
ing process? Is it then possible that an edge state might 
compete with an interior state due to CAT and/or RAT? 
Note that the four zeroth-order states in Fig. [20] are in 
fact close to one such resonance junction. The classical 
survival probabilities in Fig. [18] indicate some transport 
but further studies are required in order to character- 
ize the classical transport. The issues involved are sub- 
tle since the competition between the coexisting classical 
and quantum routes has hardly been studied up until now 
and warrants further attention. Few of the studies that 
have been performed [234L [235L 12361 . [23^ to date suggest 
that quantum effects tend to localize the along-resonancc 
diffusion on the web. 



Nearly a decade ago Leitner and Wolynes analysed [215 
the Hamiltonian in Eq. ([3]) to understand the role of the 
vibrational superexchange at low energies. Within the 
state space approach it was argued that the superex- 
change mechanisms contribute significantly to the energy 
flow from the edge states to the interior of the state space. 
More importantly the superexchange contributions mod- 
ify the IVR threshold as quantified by the transition pa- 
rameter T(E) defined in Eq. ([TO]) . From the phase space 
perspective as in Fig. [23] the direct and superexchange 
processes correspond to classical and quantum dynamical 
tunneling mechanisms respectively which are determined 
by the topology of the web. In contrast, note that the 
present studies on the model Hamiltonian in Eq. ([62]) ac- 
tually shows dynamical tunneling at fairly high energies 
and between interior states in the state space. Further 
studies on such model systems would lead to better in- 
sights into the nature of T(E). It would certainly be 
interesting to see if the criteria T(E) « I of Leitner and 
Wolynes is related to the criteria for efficient CAT and 
RAT in the phase space. Needless to say, such studies 
combined with the idea of local control demonstrated in 
the pevious section can provide important insights to- 
wards controlling IVR. The effective h is expected to 
play a central role in deciding the importance of various 
phase space structures. However intutive notions based 
entirely on this 'geometrical' criteria can be mis lead ing. 
For instance, Gong and Brumer in their studies |l46| on 
the modified kicked rotor found the quantum dynam- 
ics to be strongly influenced by regular islands in the 
phase space whose areas were at least an order of magni- 
tude smaller than the effective Planck constant. Another 



matically driven hydrogen atoms is mediated by nonlin- 
ear common resonances [23^ with phase space areas of 
the order of the Planck const ant. Furt her examples can 
be found in the recent studies I86l . 206j . 

Although for large molecules the vibrational superex- 
change viewpoint is perhaps the only choice for compu- 
tations, for a mechanistic understanding the phase space 
perspective is inevitable. Furthermore, as noted earlier, 
there exists ample evidence for the notion that even in 
large molecules the actual effective dimensionality of the 
IVR manifold in the state space is far less than the the- 
oretical maximum of (3 N — 6) . Hence detailed classical- 
quantum correspondence studies of systems with three 
or four degrees of freedom can provide useful insights 
into the IVR processes in large molecules. At the same 
time a proper understanding of and including the effect of 
partial barriers into the theory of dynamical tunneling is 
necessary for further progress. Some hints to the delicate 
competition between quantum and classical transport 
through cantori can be found in the fairly detailed studies 
by Maitra and H eller I67l | on the dynamics generated by 
the whisker map [2071 ]. Such studies are required to gain 
insights into the method of controlling quantum dynam- 
ics via su itabl e mo difications of the local phase space 
struct ures [240l 241 1. Extending the local phase space 



control idea to systems with higher degrees of freedom 
poses a difficult challenge. If the sketch in Fig. [23] is rea- 
sonable then it suggests that the real origins of the hierar- 
chical nature of IVR, due to both quantum and classical 
routes, is intricately tied to the geometry of the Arnol'd 
web. Significant efforts are needed to characterize the 
resonance web dynamically in order to have a clear un- 
derstanding of how the specific features of th e web in flu- 
ence the IVR. Recent works I6C 



61, 62, 63. 67, 68, 2201 arc 



example comes from the work by Sirko and Koch 238 
wherein it was observed that the ionization of bichro- 



beginning to focus on extracting such details in systems 
with three or more degrees of freedom which promises 
to provide valuable insights into the classical and quan- 
tum mechanisms of IVR in polyatomic molecules. Per- 
haps such an understanding would allow one to locally 
perturb specific regions of the resonance web and hence 
ultimately lead to control over the classical and quantum 
dynamics. 
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